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Abstract 

Let ^ denote space-time white noise, and consider the following stochas¬ 
tic partial differential equations: (i) ii = -|- u^, started identically at 

one; and (ii) Z = ^Z" + started identically at zero. It is well known 
that the solution to (i) is intermittent, whereas the solution to (ii) is not. 

And the two equations are known to be in different universality classes. 

We prove that the tall peaks of both systems are multifractals in a 
natural large-scale sense. Some of this work is extended to also establish 
the multifractal behavior of the peaks of stochastic PDEs on R+ x 
with d ^ 2. G. Lawler has asked us if intermittency is the same as 
multifractality. The present work gives a negative answer to this question. 

As a byproduct of our methods, we prove also that the peaks of the 
Brownian motion form a large-scale monofractal, whereas the peaks of the 
Ornstein-Uhlenbeck process on R are multifractal. 

Throughout, we make extensive use of the macroscopic fractal the¬ 
ory of M.T. Barlow and S.J. Taylor [3,4]. We expand on aspects of the 
Barlow-Taylor theory, as well. 

Keywords: Intermittency, multifractality, macroscopic/large-scale Haus- 
dorff dimension, stochastic partial differential equations. 

AMS 2010 subject classification: Primary. 60H15; Secondary. 35R60, 
60K37. 

1 Introduction, and main results 

The principle aim of this article is to answer the following question that was 
posed to us by Gregory Lawler in January of 2012 [private communication]: 

‘Research supported in part by NSF grant DMS-1307470. 
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Is “intermittency” the same property as “multifractality”? We will argue below 
that the short answer is no. Among other things, it follows that the macroscopic 
analysis of disordered systems can unravel a great deal more complexity than 
its microscopic counterparts (as compared with the theory of Paladin et al [34], 
for example). 

The two quoted terms, “intermittency” and “multifractality,” appear also in 
the title of this article. They are meant to be understood as informal descriptions 
of behavior that is commonly observed in a vast array of complex scientific 
problems in which there are infinitely-many natural length scales. This sort of 
behavior is common in, but not limited to, problems in full-blown turbulence. 

Intermittency is a well-defined mathematical property which we recall first. 
Multifractality will be treated afterward, and will require more effort. In fact, 
one of the novel parts of this paper is to set up a mathematical framework 
within which we can understand macroscopic multifractality in a way that is 
meaningful in the present context. 

In order to motivate some of the results of this paper let us consider the 
eigenvalue problem for the random heat operator, 

(i-i) 

acting on space-time functions / = f{t,x) such that f{0,x) = 1 for all a; G R, 
say. Here, M denotes the random multiplication operator, defined via 

{Mf){t,x) := f{t,x)^{t,x) (t>0,xGR), (1.2) 

where ^ is a space-time white noise. Because ^ is a random Schwartz-type 
distribution, and not a nice classically-defined random process, the operator M 
needs to be understood in integrated form: As a Wiener-integral map when / 
is nonrandom; and more generally as a Walsh-type stochastic integral operator 
when / is a predictable random field. 

It is well known, and also easy to see, that the spectrum of H is all of 
R, and that the eigenfunction of H that corresponds to eigenvalue A G R 
is ex.p{—Xt)u{t, x), where u = u{t,x) solves the parabolic Anderson model, 
which we may understand rigorously as the solution to the following Walsh- 
type stochastic partial differential equation [39]: 

u=hu" + ul on(0,oo)xR, 

2 , , (PAM) 

subject to m( 0 , a;) = 1 for all a; G R. 

It is well known that the parabolic Anderson model (PAM) has a unique 
solution which is a predictable random field that is continuous in both variables 
t and x} Moreover, there exist finite and positive constants Uo,U, Lg, L such 
that 

Lq exp < E (|M(t , a;)!^) ^ Uq exp (Uk^t) , (1.3) 

^ The term predictable is used as in Walsh [39], and refers to predictability with respect to 
the filtration generated by the white noise 
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valid uniformly for all real numbers t > 0, k ^ 2, and x G R. [In fact, these 
moments do not depend on the value of x.] See for example Chapter 5 of [28] for 
a self-contained account. As a consequence, the moment Lyapunov exponents of 
the solution u are strictly positive and finite, where the kth moment Lyapunov 
exponent of u is defined as 

X{k) := lim logE (|■u(^ , 0)1*) , (1.4) 

t—^OO ^ ' 

for every real number fc > 0. 

One can devise a subadditivity argument in order to show that these Lya¬ 
punov exponents exist. It is widely believed that the following Kardar formula 
holds [25-27] :2 

X{k) = — l)/24 for all real numbers k > 0. (1-5) 

The a priori bounds (1.3) and convexity considerations together imply that 

k !->■ k~^X{k) is strictly increasing on [2 ,oo). (1.6) 


This property is known as intermitteney and was referred to earlier on in a 
broader context. 

It is possible to argue that because of intermitteney the solution to the 
parabolic Anderson model (PAM) has very tall peaks on infinitely-many differ¬ 
ent length scales. We have not yet described very precisely what it means to 
“have tall peaks on infinitely-many different scales.” Still, the adjective “mul¬ 
tifractal” is supposed to reflect the presence of such behavior. The interested 
reader can find different heuristic accounts of how intermitteney might explain 
multifractality in the Introductions of Bertini and Cancrini [5] and Carmona and 
Molchanov [7]. Chapter 7 of Khoshnevisan [28, §7.1] yields a related but slightly 
more precise explanation. A recent paper by Gibbon and Titi [20] contains an 
account of how, and why, intermitteney and other attempts at describing “mul¬ 
tifractality” arise naturally in a large number of multi-scale problems in science. 

Now that we recalled “intermitteney,” within context, we begin to propose 
a mathematical model of “multifractality” that is general enough for our later 
needs, but concrete enough that it is amenable to exact analysis. It turns 
out to be helpful to not begin with a too-concrete setting such as stochastic 
PDEs. Therefore, let X := {A(t)}tgRn be a real-valued stochastic process with 
continuous trajectories. We are interested in saying that X has tall peaks that 
are “multifractal” in a large-scale sense. In order to do this, let us posit that 
the tall peaks of the stochastic process X are described by a gauge funetion 
g : R+ —>■ R+. By this we mean that (/ is a nonrandom and increasing function 
such that limr^oo5('c) = oo and 


lim sup 

||t||->oo 


\X{t)\ 

5(PII) 


= 1 


a.s. 


(1.7) 


^Khoshnevisan [28, Chapter 6] has shown that L ^ — l)/24 for every integer fc ^ 2. 

Thus, in particular, A(fc) ^ k{k‘^ — l)/24 for all integers k ^ 2. Borodin and Corwin [6] have 
found the first completely-rigorous proof of (1.5) in the closely-related case that /c ^ 2 is an 
integer and the initial value of (PAM) is replaced by 5o- 
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That is, we assume that there exists a non-random gauge function g which 
describes the largest possible scale on which we can compute the tall peaks of 
X in a macroscopic sense. Let us fix the gauge function g in our mind, and then 
consider the random set, 

Vx{l) ■■= e : \\t\\ > CO, M , (1.8) 

where 7 > 0 is a tuning parameter and Cq is a fixed positive constant. We 
may think of every point t G Vxil) as a tall peak of the process X —suitably 
normalized by g — viewed in length scale 7 . 

Because Vxih) Q 'Px(l 2 ) when 71 > 72 , we might expect the existence 
of a phase transition that does not depend on our prior choice of cq. This is 
correct, and not difficult to explain. In fact, the growth condition (1.7) implies 
that there exists a uniquely-defined transition point at 7 = 1: Vxil) is a.s. 
unbounded when 7 < 1 whereas Vx{l) is a.s. bounded when 7 > 1 . Motivated 
by this simple observation, let us introduce the following. 

Definition 1.1. Choose and fix a gauge function g and consider the tall peaks 
of X, as defined in (1.8). We say that the tall peaks of X are multifractal 
when there exist infinitely-many length scales 71 > 72 > • • • > 0 such that with 
probability one, 

Dimj, {Vx{li+i)) < Dimjj {Vx{li)) for all i > 1 . (1.9) 

If, on the other hand. 


DimH(T’x(7)) = 

then we say that the tall peaks of X are monofractal? 

In the preceding, Dim^ [E) denotes the macroscopic Hausdorff dimension of 
E C R". This notion is due to Barlow and Taylor [3,4], and will be reviewed 
in the next section. For the time being, it suffices to know only that Dim^ {E) 
is a real number between 0 and n, and describes the large-scale geometry of 
E C R” in a way that parallels how the usual notions of fractal dimension try 
to describe small-scale geometry. 

In order to see how our definition of multifractality can be used in stochastic 
PDFs, let us return to the stochastic heat equation, (PAM). A deep theorem of 
Mueller [30] asserts that 

u{t ,x) > 0 for alH ^ 0 and x G R, (l-H) 

®In principle, it can happen that the tall peaks of X are neither multifractal nor monofrac- 
tal. That is, when (1.9) holds for finitely-many 71 > 72 > ■ • • > 7 n- It might be interesting 
to construct non-trivial examples of stochastic processes whose tallest peaks are of this latter 
type. We aren’t aware of any natural examples at this time. 


whenever 7 < 1 , 
whenever 7 > 1 , 


(I.IO) 
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outside a single event of probability zero. Because the natural logarithm is 
strictly monotonic, the tall peaks of the random function x i—>■ u(t, x) are in 
one-to-one correspondence with the tall peaks of a; i—>■ Xt{x) at all times t > 0, 
where 

/32\ 

:= logu{t,x). ( 1 . 12 ) 

Conus et al [12] have proved that for all t > 0, 


^ Xt(x) 

0 < lim sup — 
x-j-oo g{x) 


< oo 


a.s., where g(x) := (log^x)^/^. 


(1.13) 


More recently, Chen [8] has found the following improvement to (1.13): For all 
t > 0 , 

lim sup ^ = 1 a.s. (1-14) 

x^oc g(x) 

In other words, g(x) := (log^ x)^^^ is a natural gauge function for measuring 
the tall peaks of the random height function x i—>■ Xt{x). The particular normal¬ 
ization of X in (1.12) is there to merely ensure that the lim sup is one, which 
matches the form of (1.14) with that of (1.7). Having said this, we are ready 
to mention one of the results of this work. 

Theorem 1.2. Consider the gauge function g{x) := (log^ x)^/^ and, for all 
t > 0, define Xt = {Xt{x), x € R} as in (1.12). Then the tall peaks of Xt are 
multifractal with probability one for all t > 0. In fact, for every t,j>0, 

Dim„ (^^^,( 7 )) = 1 - as-, (1-15) 


where Dim^ (E) < 0 means that E is bounded. 

Theorem 1.2 shows that the tall peaks of the solution u to (PAM) are multi¬ 
fractal in the sense described more precisely in Theorem 1.2. And we can think 
of (1.15) as a description of large-scale “multifractal spectrum” of the peaks 
of the random field u. Theorem 1.2 is a consequence of a much deeper result 
which says, using informal language, that the tall peaks of the solution to every 
known multifractal nonlinear stochastic PDF [18] are generically multifractal. 
See Theorem 7.2 and Corollary 7.3 below. 

Earlier we mentioned that u is intermittent in the sense of (1.6). Together 
with Theorem 1.2, this assertion says that intermittency and multifractality can 
coexist. Next we describe a non-intermittent system whose tall peaks are also 
multifractal. Consequently, non-intermittency and multifractality can also co¬ 
exist, whence follows the negative answer to Professor Lawler’s question with 
which we began. 

Consider the linear stochastic heat equation, 

Z=^Z'' + f on(0,oo)xR, 
subject to Z(0 ,x) =0 for all x G R. 


5 



It is well known that a solution exists, is unique, and is a continuous centered 
Gaussian random field. It is also possible to prove that the moment Lyapunov 
exponents of Z are all zero, and so Z is not an intermittent process."^ 

Also, one can prove fairly easily that for all t > 0 the following holds with 
probability one: 

limsup = (t/7r)^/^ a.s., where g(a;) := (2 log , (1-16) 

x^oo g{x) 

A “steady state” version of this fact appears in print, for example, in Collela 
and Lanford [10, Theorem 1.1(c)]. And the fact itself follows by specializing an 
even earlier, very general, theory of Pickands [36]; see also Qualls and Watanabe 
[37]. 

Let 

Xt{x) := ,x). (1.17) 

The following shows that the tall peaks of Xt —hence also those of Z{t, •)—are 
multifractal even though Z is non-intermittent. 

Theorem 1.3. Consider the gauge function g{x) := (log^a;)^/^ and, for all 
t > 0, define Xt = {Xt{x), x € R} as in (1.17). Then the tall peaks of Xt are 
multifractal with probability one for all t > 0. In fact, for every t,j>0, 

Dim^ (^^^,( 7 )) = 1 - 7 ^ a.S; (1.18) 

where Dim^ (E) < 0 means that E is bounded. 

Theorems 1.2 and 1.3 are both consequences of two more general theorems 
about multifractal random fields [Theorems 4.1 and 4.7 below]. Those general 
theorems have other interesting consequences as well. Let us mention one such 
result. 

Theorem 1.4. Let B denote a one-dimensional Brownian motion and U{t) 
:= exp(—t/2)i3(e*) an Ornstein-Uhlenbeck process on R. Then, the tall peaks 
of B are a monofractal, whereas those of U are multifractal in the following 
sense: For every 7 > 0, 


Dimjj 


s > e® 


B{s) 

( 2 s log logs)i /2 



1 */7<l, 
0 i/7>l, 


(1.19) 


almost surely, whereas 

' (2 log‘.72 > 7} = 1 - (1.20) 

almost surely, where we recall Dimjj (E) < 0 means that E is bounded. 

^This can be deduced from inspecting the proof of Theorem 2.1 of ref. [18], for example. 
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Theorem 1.4 will be proved in two parts: (1-19) is proved below in Theorem 
3.2; and (1.20) is proved in Theorem 5.1. 

We have included a final §8 wherein we state and prove a suitable variation 
of Theorem 1.2, in which the stochastic partial differential equation (PAM) is 
replaced by a similar-looking object with u" replaced by Au, the Laplacian, the 
latter acting on a space variable x G IR"^ in place of R. And space-time white 
noise is replaced by a centered Gaussian noise that is white in time and suitably 
correlated in space to ensure the existence and uniqueness of a well-tempered 
solution. Among other things, such equations are well-known models that play 
a role in the large-scale structure of the universe; see §8 for more details. 

Let us conclude the Introduction with a word on notation. From now on 
we will write (pt (x) and (p(t , x) interchangeably, depending on which is more 
convenient, for any space-time function [or even generalized function/noise] </>. 
This choice ought not inconvenience the reader, since pt is never used to denote 
the time derivative of </>; rather, it is the standard probabilistic notation for 
describing the time evolution of </>. 

2 Dimension and density 

Let us begin by recalling the Barlow-Taylor definition of the macroscopic di¬ 
mension [3,4] Dimjj A of a set E C R"^. 

Define, for all integers n ^ 0, 

V„ := [-e" , e")”^, 5o := Vo, and 5„+i := V„+i \ V„. ( 2 . 1 ) 

We might sometimes refer to as the nth shell in R'^. 

An important idea of Barlow-Taylor [3,4]—see also Naudts [32] for a pre¬ 
cursor to this idea—is to construct a family of Hausdorff-type contents on each 
shell, and then use the totality of those contents in order to define a family of 
Hausdorff-type contents on all of R'^. Once this is done, a notion of macroscopic 
Hausdorff dimension presents itself quite naturally. 

Definition 2.1. Let B denote the collection of all sets of the form 

Q{x ,r) :=[xi,xi+r) X ■■■ X [xd ,Xd + r), (2.2) 

as X := (xi ,..., Xd) ranges in R'^ and r in (0 , oo). If Q := Q{x , r) is an element 
of R, then we may refer to Q as an upright box with southwest corner x and 
sidelength side(( 5 ) := r. 

Choose and fix some number cq > 0, and define for every set E C R"^, all 
real numbers p > 0 , and each integer n ^ 0 , 

(2.3) 

where the infimum is taken over all upright boxes Qi,..., Qm of side ^ cq that 
cover A n . We may think of t')) {E) as the restriction to the nth shell of 
the scaled p-dimensional Hausdorff content of E. 
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Definition 2.2. The Barlow-Taylor macroscopic Hausdorff dimension oi E 
is 


Dinijj E := Dimjj(_E) := inf > 0 : ^p{E) < oo j . (2-4) 

By (2.4), any bounded set E C R'^ has Dinijj E = 0. We will leave the 
following simple fact as exercise for the interested reader. 


Lemma 2.3. The numerical value o/Dim^ E does not depend on cq > 0. 

Thus, we could choose cq = 2 or cq = e/-\/2, in our definition of t'", without 
affecting the value of Dim^ E. It is important to point out though that cq = 0 
can lead to a different value of Dim^ E. 

In order to see why we are ruling out the possibility of cq = 0, let us define 
for all E C R"^ and r > 0 a new set E^{r) C lA as follows: 

E^{r)-.= {x&'il'^-. EnQ{x,r)^0). (2.5) 


The notation is basically due to Barlow and Taylor, who observed the following 
[3, Lemma 6.1] but stated it using slightly different language: For all E C R'^ 
and r > 0, 

Dimjj Fi = Dimjj(£i^(r)), (2.6) 

where Dim^ on the right-hand side of (2.6) is the discrete Hausdorff dimension 
on of Barlow and Taylor [3,4]. In other words, because we chose only covers 
of E that have side > co > 0, the local structure of E does not affect the value 
of its macroscopic Hausdorff dimension. Put yet in another way, this shows that 
the Barlow-Taylor definition of Dim^ E quantifies the large-scale geometry of 
E, without obstruction by the microscopic structure of the set E. 

Since the particular value of cq > 0 does not matter, from now on we follow 
the choice of Barlow and Taylor [3,4], and set 


Co = 1. 


(2.7) 


The following is a macroscopic counterpart of a familiar result about micro¬ 
scopic Hausdorff dimension. 

Lemma 2.4. Suppose f : E ^ R^* is a Lipschitz function on E G R*^ and 
satisfies the growth condition liminf |/(a;)|/|a;| > 0. Then, 

x^E, |ai|— 

Dimjj f{E) ^ Dimjj E. (2-8) 

In particular, if f : R'^ —>■ R*^ is bi-Lipshitz on E C R*^; that is, there exists a 
positive constant L ^ 1 such that 

L~^\x - 2/K l/(a;) - f{y)\ < L\x - y\ for all x,y € E. 


Then 


Dimjj f{E) = Dimjj E. 


(2.9) 



Proof. The proof is similar to that of the same assertion for ordinary Hausdorff 
dimension: For every p > Dim^ E, we can find upright boxes Qj^n G S —indexed 
by 1 ^ ^ mn, n ^ 1 —with side{Qj^n) ^ 1 such that: 

(i) if n C for all n ^ 1; and 

(ii) EZiETM^^MQj.n)/e-y <oo. 

It is clear that f{E) C fiQj,n)- Since / is Lipschitz on E we can 

hnd a finite and positive constant q such that 


\f{x)- f{y)\^q\x-y\ for all x,y G S. (2.10) 


Without loss of generality, we may assume that g ^ 1 is an integer; otherwise, we 
replace g by 1 + [gj everywhere. In particular, it follows that every f{Qj^n) can 
be covered with an upright box whose sidelength is between 1 and qside{Qj^n)- 
In addition, the growth condition on / implies that there is a constant £ > 0 
such that f{x) ^ £e” for every x £ Eld Sn This implies that 


OO 

E 

n=0 




n=0 


j=i 


/ side(Qyn) 


ee" 


for all n ^ 1, 


( 2 . 11 ) 


whence < oo by (ii). This proves that Dimjj f{E) < p for all 

p > Dimjj E and implies (2.8). 

Finally, the bi-Lipshitz condition implies that both / and its inverse f~^ on 
f{E) satisfy the growth condition. Hence, we make two appeals to the first part 
of Lemma 2.4, once for / and once for /“^, to see that (2.9) holds. □ 

Lemma 2.4 is the large-scale/macroscopic analogue of the following well- 
known fact: If / : if —>■ R is locally Lipschitz continuous, then 


dimg f{E) < dimjj E, (2.12) 

where “dim^” denotes the usual [microscopic] Hausdorff dimension in R*^. Let 
us, however, observe that (2.8) does not hold when / is only Lipschitz. For 
example, set f{x) := ln(x) for a; ^ 1 and E = exp(lN) to see that 


Dimjj /(if) = Dimjj M = 1 > 0 = Dim^ (exp(M)). (2-13) 


In the above, IN := {1,2 ,...} denotes as usual the set of all natural numbers. 
We will present an interesting example of E in Remark 5.2 below which shows 
that ( 2 . 8 ) does hold for f{x) := ln(a:) even though / does not satisfies the growth 
condition in Lemma 2.4. 

Next, let us mention a technical estimate, which is a “density theorem.” The 
following is a large-scale analogue of the classical Frostman lemma, and basically 
rephrases Theorem 4.1(a) of Barlow and Taylor [3] in a slightly different form 
that is more convenient for us. 
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Lemma 2.5 (A Frostman-type lemma). Choose and fix an integer n ^ 1, and 
suppose E C Sn is a Borel set in Let p, denote a finite non-zero Borel 
measure on E, and define for all p 0, 

Kp := : Q e B, Q C Sn, side(Q) ^ l| . (2.14) 

Then, v^{E) > p{E). 

Remark 2.6. The constant Kp of Lemma 2.5 typically depends on n as well, 
and is always finite and positive. 

Let be a Borel measure on IR,'^. For all E C IR"^ let Den,y E denotes the 
upper density of E with respect to v. That is, 


Deiii^ E := Deny(i?) := limsup 

t—¥00 


v{EC[-t,tY) 

w 


(2.15) 


The following describes an easy-to-verify sufficient condition for a set E in 
to have full macroscopic Hausdorff dimension. 

Lemma 2.7. Let v denote either the Lebesgue measure on R'^ or counting 
measure on a sublattice e'lfi of R*^ for some e > 0. Then, for all Borel sets 
E C R'^, if Denj/ E > 0, then Dim^ E = d. 

Proof. We will prove the lemma only in the case that iz denotes the Lebesgue 
measure. When iz is the counting measure on the result is proved in almost 
exactly the same way. 

Barlow and Taylor [3,4] have introduced another large-scale notion of Haus¬ 
dorff dimension, which we write as follows: 


Dim E := Dim (E) := inf > 0 : lim v’^iE) = o} • (2-16) 

^ ^ L n^oo ^ J 


It is easy to verify that Dim ^^ E ^ Dim^ E for all E C R'^; therefore we might 
think of Dinij^ as the lower Hausdorff dimension, in the macroscopic sense. Our 
goal is to prove the following somewhat stronger statement: 


Denj/ A > 0 Dim ^^ E d. (2-17) 

It follows from (2.17) that Dim^ E ^ d. Since Dimjj E ^ Dimjj(R‘^) = d (see 
Barlow and Taylor [4, Example 4.1]), (2.17) completes the proof. 

For every real number a > 1 and integers n 0 dehne 

Vn(a) := [-a",a”)^ 

5o(a) :=Vo(a), (2.18) 

5„+i(a) := V„+i(a) \ V„(a). 
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Furthermore, 

Ai(K(a)) ^ |V„(o)| = 2^0”“^ for all n ^ 0, (2.23) 

and 

M(V„+i(a)) ^ 2^(Den, E - (2.24) 

for infinitely-many integers n ^ 1, where {(5n}))^i is a sequence that satisfies 
lim„^oo <^n = 0. Therefore, 

^Ji{Sn+l{a)) = ^l{Vn+l{a)) - KVn{a)) ^ (a^Den,.E-(5„-l)2V^ (2.25) 

for infinitely-many integers n ^ 0. It is easy to adapt Lemma 2.5 to an analogous 
statement about i^p{E; a) for all choices of a [not only a = e]. That endeavor 
shows us that i^^{E; a) ^ a‘^ Deui, if — — 1 for infinitely-many integers n ^ 0. 

Therefore, 

limsup a) > 0, (2.26) 

n—¥oo 

provided that we choose a > max{l, (Den,^ This and (2.20) together 

imply (2.17) and hence the lemma. □ 


11 



3 Peaks of Brownian motion 


Consider the set of times, after time t = exp(e) [say], at which the Brownian 
motion has LIL-type behavior. That is, let us fix some parameter 7 > 0, and 
consider the random set 



B{s) 

( 2 s log logs)^/^ 



(3.1) 


We are using a notation that is consistent with that in (1.7) and (1.8), where 
g{x) := ( 2 a;log_,_ log_,_ a;)^/^ is the gauge function that comes to us naturally 
from the standard law of the iterated logarithm for Brownian motion. That is, 
the following: 

Proposition 3.1 (The law of the iterated logarithm). With probability one: 

1- 'Pb{i) is unbounded a.s. when 7 ^ 1 ; 

2. 'Psil) is bounded a.s. when 7 > 1 . 

Actually, the standard textbook form of the LIT refers only to the case that 
7 < 1 and 7 > 1. The critical case 7 = 1 follows from Motoo’s work [29, 
Example 2]. The following shows that when Vb{i) is unbounded, that is when 
7 ^ 1 , it is a macroscopic fractal of dimension one. [Of course, Vb{i) does not 
have a remarkable macroscopic structure when 7 > 1 .] 

Theorem 3.2. Assertion (1.19) of Theorem I .4 holds; that is, 


Dimjj Pb(7) = 1 a-s. for all ^ 1. 


(3.2) 


In the subcritical case where 7 < 1, Theorem 3.2 follows immediately from 
Lemma 2.7 and the next statement. 

Proposition 3.3 (Strassen [38]). Let v denote the Lebesgue measure on 
Then the following assertions are valid a.s.: 

1. Den,y 'Pb{i) > 0 when 7 < 1 ; and 

2 . Den,y ’Rb(I) = 0. 

Proof. One can easily adapt a result of Strassen [38, eq. (11)] about random 
walks to a statement about linear Brownian motion in order to see that with 
probability one. 


Den,,PB(7) = 1 - exp 





(3.3) 


as long as 7 ^ 1. This does the job. 


□ 
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Proposition 3.3 shows that Theorem 3.2 is interesting mostly in the critical 
case. In the critical case, the random set Pb( 1) is comprised of tall peaks of 
maximum possible asymptotic height. And Theorem 3.2 shows that the set 
of tall peaks of critical height has full dimension although it has zero density 
[Proposition 3.3]. 

With the preceding remarks in mind, let us consider the following random 
Borel measure that is supported in Pb(1): 


^(G) := |Pb( 1) n G| for all Borel sets G C [4 , oo). 


(3.4) 


where j • • • j denotes the 1-dimensional Lebesgue measure. The following is the 
key step in the proof of Theorem 3.2. 

Proposition 3.4. ^ 

We will begin our proof of Proposition 3.4 shortly. But first, let us apply 
this proposition in order to establish Theorem 3.2. 

Proof of Theorem 3.2. The preceding remarks tell us that we need to only con¬ 
sider the critical case, 7 = 1 . Because ^jl[x,x -|- r) ^ r, it follows that Ki = 1, 
where Kp was defined in (2.14). Lemma 2.5 and Proposition 3.4 together imply 
that 


^ bl‘5„) ^ ^e ”/x(5„) = 00 a.s. (3.5) 

n=4 n=4 

In particular, it follows that Dim^ P_b( 1) ^ 1 a.s., which is the desired result. □ 
In order to derive Proposition 3.4, let us consider the events 

£t := ^uj G Tl : B{t){u}) > (2tloglogt)^/^| for all t > 4. (3.6) 

It is easy to see from I’Hopital’s rule that if X has the standard normal distri¬ 
bution, then P{X > z} is to within a multiplicative constant of z~^ exp(—z^/ 2 ) 
uniformly for all z ^ 1. The following is a consequence of this fact and the strict 
positivity and the continuity of the Gaussian density function: There exists a 
finite constant c > 1 such that 


1 


clogt(log logt)^/^ 


< p (St) ^ 


log<(log log t)l/^ 


for all t ^ 4. 


(3.7) 


Next we estimate P(£lt | £s) for various choices of s < t. The first quantifies 
the well-known qualitative fact that £s and £t are approximately independent 
when t ^ s. 


Lemma 3.5. For all real numbers t > s > A, 

t > 4s(loglogs)(loglogt) ^ P(£lt |£ls) ^ cP(£lt), 
where c G (0, 00 ) does not depend on {s ,t). 


(3.8) 
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Proof. We recall the following well-known bound, which is essentially Lemma 
1.5 of Orey and Pruitt [33]: If U and V are jointly distributed as a bivariate 
normal with common mean zero, common variance one, and covariance p, then 
there exists a finite constant c such that 

P([/ > a\V > 5) < cP{C/ > a} whenever \p\ < {ab)~^. (3.9) 

Next we apply the preceding by setting U := V := a := 

(2 log log and b := (2 log log Note that p = satisfies 0 < p < 

{ab)~^ because t > 4s(loglogs)(loglogt). The lemma follows from (3.9). □ 

When t and s are not too far apart, we do not expect P{£t | iPs) to have the 
same order of magnitude as P (£!*). The following provides us with a quantitive 
estimate of P{£t | i^^s) in this case. 

Lemma 3.6. There exists a finite eonstant c such that for every t > s > 4, 

^^^-l^-)^ (logs)('--)/(4W (3.10) 

Proof. According to Lemma 1.6 of Orey and Pruitt [33], if U and V are jointly 
distributed as a bivariate normal with common mean zero, common variance 
one, and covariance p, then there exists a finite constant c such that 

P(C/ > a\V > a) ^ cexp ~ a ^ 0. (3.11) 

We apply this inequality with U := t~^^‘^B(t), V := s“^/^i?(s), and a := 
(2 log log s)^/^ to find that 

Fi£t\£s) < P (B(t) > (2t log logs) I £:,) = P([/ > a I P > a). (3.12) 

Thus, the lemma follows from (3.11). □ 

We are prepared to verify Proposition 3.4. 

Proof of Proposition 3.4. For all iV ^ 4 define 

N 

^^:=^e-X5„). (3.13) 

n—4 

We intend to prove that Soo ■= lim^v-^oo S'at is infinite a.s. Because ES'at = 

J2n=4 /e”-i ds, we may apply (3.7) in order to see that 

ES'jv ^ Co(logA^)i/2 for all N 4, (3.14) 

where cq is a positive constant that does not depend on N. Next, we estimate 
the variance of S'at. Clearly, 

E(5^)^Qi+Q2, (3.15) 
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where 


N 


n—4 


Qi:=2j2 e"^” 11 ]P(fs n St) ds dt, 

/en —1 <5<^<en 


n m 

Q 2 :=2 EE e-”-™ r ds / dt P(i?s n St). 


(3.16) 


4^n<m^N 

The elementary bound V{SsC\St) ^ P(iPs) yields 


Qi < 2E5Ar. 


(3.17) 


We estimate Q 2 by splitting the double sum according to how much the 
summation variable m is greater than the summation variable n. Before we 
hash out the details, let us first note that, according to Lemma 3.6, whenever 
m > n ^ 4, 


f ds f dt P(£ls n St) 

^ Cl ^ ds dt P(£ls) exp log log s^ 

^ Cl f ds f dt P(£ls) exp (—02 

Je”-! \ 


t - e" 


(3.18) 


^00 m 

^cie™/ ds dt P(fs)e-"^‘^°s” = ^ / P(f«)ds, 
Jen-1 Jo lognJe—i 


where ci, 02,03 are finite and positive constants that do not depend on {n,m). 
Consequently, 


EE e 

4^n<m^(n+Q; log n)AN 


—n—m 


f ds f dtP(£lsniPt) 

Jen-1 Je™-! 



(3.19) 


cocxUSjn , 


for all a > 0. We emphasize that ci, 02 , and 03 are finite constants that do not 
depend on (N , a), and a is, so far, an arbitrary parameter. 

Next we observe that there exists a > 0 large enough such that for all 
integers m,n ^ 4, 

m>n + alogn e™"^ > 4e"(logn)(logm). (3.20) 


We will choose a to be this particular value, both in the preceding and in what 
follows. In this case, it is then easy to see that t > 4s(loglogs)(loglogt) for 
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every s G and t € Sm, as long as n, m ^ 4 are arbitrary integers that satisfy 
m > n + alogn. Thus, Lemma 3.5 ensures that, for this particular choice of a, 


EE 


4^n<n+a log n<m^N 


ds 


dt P(fs n£’t) 



(3.21) 

C4 |ES'jv|^ , 


where C 4 is a finite constant that does not depend on n. Thanks to (3.16), 
(3.19), and (3.21), it follows that Q 2 ^ csoES'at + C 4 |E 5 'Arp, uniformly in all 
N ^ 4. Thus, it follows from (3.14) and (3.17) that 

E(5^) =0(|E5^,|^) as N ^ 00 . (3.22) 

Since ES'^r t 00 as —>■ 00 [see (3.14)] and Soo ^ Sn for all N, 

E {-S'oo = 00 } ^ liminf P {Sn ^ ^ES'at} ^ i liminf > 0, (3.23) 

N^oo '■ ^ J 4 N^oo E(S^^) 

thanks to the Paley-Zygmund inequality [35, Lemma 7 ], which states that if Z 
is in L^(P) and Ill'll 00 > 0 , then 

P{Z>lEZ)si||ll. (3.24) 

An appeal to the Hewitt-Savage 0-1 law completes the proof. □ 


4 General bounds 


Let X := {Xt}t£T be a real-valued random field with continuous trajectories, 
where T C is either one of the 2^^ standard closed orthants of R*^, or T is R'^ 
itself. 

For all real numbers & G (0, 00 ) we can define 



c(6) := — limsupi; ^suplogP {Xt > 2 ;} , 

z->oo teT 

(4.1) 

and 


C{b) := — lim inf z ^ inf log P {Xt > z} . 

2 —>-oo tGT 

(4.2) 

Of course. 

0 ^ c(b) ^ C(6) ^ 00 for all 5 > 0. 


Define 

t to be the £°°-norm of t G R*^; that is. 



t := max \ti\ for all t := {ti,... ,td) € R*^. 

l^i^d 

(4.3) 
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Then it should be intuitively clear that, under mild conditions on large-scale 
smoothness and asymptotic pairwise independence of X, one ought to be able 
to prove that 


0 < limsup 

|t|—>-oo 


(log|t|)i/f> 


< 00 


a.s., 


(4.4) 


provided additionally that 0 < c{b) < C{h) < oo for a certain special value of 
5 G (0, oo). In other words, we might expect that if 0 < c{b) ^ C{b) < oo, then 
the tall peaks of the process X are typically gauged, to within a constant, by 
the function t >->■ (log 

The main results of this section are two general macroscopic Hausdorff di¬ 
mension estimates; see Theorems 4.1 and 4.7 below. The first theorem describes 
conditions, similar to those outlined earlier, which ensure the upper bound in 
(4.4), and also bound from above the macroscopic Hausdorff dimension of the 
set of times that Xt exceeds a [correct] constant multiple of (log \ The sec¬ 

ond theorem turns out to be a much more subtle result that produces matching 
lower bounds for the dimension of these exceedance times. Thus, we begin 
with a general upper bound for the Hausdorff dimension of the tall peaks of a 
stochastic process. 


Theorem 4.1 (A general upper bound). Suppose that there exists b G (0,oo) 
such that c{b) > 0 and for all 7 G (0 , d), 


wGT I 


c{b) 


i/b' 


sup P sup Xt > ( —logs ) ^ < s as s ^ CO. (4.5) 


Then, 


Furthermore, 


i/b 




(4.6) 


i/b^ 

DimH-^tGT: |t| > exp(e). At ^ ^ log |t| ) < d - 7 , (4.7) 


for all 7 G (0 , d). 

Proof. It suffices to prove the result in the case that T = 1R,^|.. The other cases, 
including T = follow from this after making small adjustments. 

The stated limsup result is a more-or-less standard exercise in the Borel- 
Cantelli lemma, and the upper bound on the dimension follows from a routine 
covering argument, the likes of which are familiar for bounding the microscopic 
Hausdorff of a random set. We include the proof for the sake of completeness. 
Let us write Dim^ E < 0 when E is a bounded set. We plan to prove that 

(4.7) holds for all 7 > 0; (4.6) follows immediately from this formulation of 

(4.7) . From now on we choose and fix an arbitrary 7 > 0. 
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Our goal is to prove that 


DirUg Au ^ c? — 7 a.s., 


(4.8) 


where 

A„ := |t G ]R,+ : \t\ > exp(e), Xt > log |t| j | . (4.9) 

Condition (4.5) ensures that for all e G ( 0 , 7 ) there exists a finite constant Kg 
such that for all m := (mi,..., nid) G RiJ. that satisfy ||m|| > exp(e), 


P{A„ n Q(m, 1) 0 } < P < sup Xt^ 


—— log \m 
c{b) 


i/b' 


< 


Kg 


\mp~ 


(4.10) 


Therefore, we can cover A„ fl Sn by upright boxes of sidelength r = 1 in order 
to see that for all p > 0 and n ^ 0 large, 


Q(m,l)c5n, 

^ Kge-'^p E 1 ^ 1 “'^+^ 

Q(m,l)c5n 


(4.11) 


Whenever n G N and m G are such that Q{m ,1) lies entirely in 5„, then 
it must be that |m| ^ exp(n — 1). Since contains at most const • exp(n(i) 
upright boxes of sidelength one, it follows that 


E [:/”(A„)] < const • 

for all n ^ 1 sufficiently large. In particular. 


E 


00 


.n—0 


:(A„) 


<00 ifp>fi — 7 + e. 


(4.12) 


(4.13) 


This proves that Dimjj A„ ^ p a.s. for all p > d — 7 + e. Send p d — 7 + £ and 
then e 4 , 0, in this order, to deduce (4.8) and hence the theorem. □ 


We now move on to a perhaps more interesting study of lower bounds for 
Dimjj of high peaks of X] more specifically, our next result shows that one 
can sometimes obtain good lower bounds on the macroscopic dimension in the 
statement of Theorem 4.1. 

Consider, for every 7 G (0, d), the random set 

\ l/b 
logt) 



Ai:= it eT 


> exp(e), Xt ^ 


Cib) 
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We plan to show that, under some conditions on the process X, 

DimjjA^^d —7 a.s., (4-15) 

thus obtaining a complimentary bound to that of Theorem 4.1. 

The standard way to obtain lower bounds on the ordinary Hausdorff dimen¬ 
sion of a set is to find a smooth measure on that set; see Lemma 2.5 for example. 
The said smooth measure is typically “uniquely canonical,” and readily guessed 
when the set in question has good local structure. The remaining work is in 
determining the exact order of the smoothness of the canonical measure. 

In the present setting of macroscopic dimension, one can in fact prove that 
the random set does not support a “uniquely canonical” measure. Instead 
we have to use a different sort of argument. 

Before we proceed we need to develop a few general facts about macroscopic 
Hausdorff dimension. We will also introduce some notation that will be used 
throughout the remainder of this section. 

Definition 4.2. Let n„ be a finite collection of points in for every integer 
n ^ 0. Given a real number 0 G (0,1), we say that {n„}5^g is a 9-skeleton of 
if there exists an integer N = N{9) such that: 

1. For every n ^ N , 

y g(a:,e®”) C5„; (4.16) 

xGlln 

2. If X and y are distinct points in n„ for some n'^ N ^ then 

Q(a;,e^")ng(y,e^”) = 0; and (4.17) 

3. There exists a constant a = a{d, 9) G (0,1) such that 

where “| • • • |” denotes cardinality. 

Given some 6* G (0,1), R*^ has uncountably-many 0-skeletons. From now on, 
we choose and fix one such choice, and denote it by n(0) := {n„(0)}))TQ. For 
instance, we might wish to opt for the following construction, to be explicit: 

n„(0) := An{9) X • • • X An{9) [d times]; (4.19) 

where 

An{9):= U {e"+je''"} (4.20) 

iez 

Other constructions are also possible. The property of n(0) that we are after 
this: n„(0) is basically a “near-optimal e^"-packing” of for all large n. Part 
2 of Definition 4.2 codifies the precise notion of “packing” and Part 3 makes 
precise our “near-optimal” sense. 
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Definition 4.3. Let E C be a set and choose and fix some real number 
0 G (0,1). We say that E is 9-thick if there exists an integer M = M(6) such 
that 

EnQix,e^'^) 0, (4.21) 

for all X G n„(6*) and n ^ M. 

In words, 0-thick sets are “approximately self-similar sets.” 

We are ready to present one of the novel technical statements of this section. 

Proposition 4.4. If E C is 9-thick for some 0 G (0,1), then 

Urnr^E'^ d{l-9). (4.22) 

Proposition 4.4 presents us with a strategy for obtaining a lower bound for 
Dimjj F that can be different from a Frostman-type method [Lemma 2.5]. The 
following is an immediate consequence of Proposition 4.4. 

Corollary 4.5. If F C R'^ has a 9-thick subset E for some 9 G (0,1), then 
Dimjj F ^ c?(l — 0). 

In principle, our definition of 0-thickness depends on our a priori choice of 
a 0-skeleton 11(0). We are not aware of any example where the choice matters 
very much. But even if the choice did matter, we can interpret Corollary 4.5 
as saying that we can always obtain a lower bound on Dim^ F by finding a 
0-skeleton 11(0) of R*^ and a set if C F that is 0-thick with respect to our choice 
of skeleton. In that case, Dim^ F ^ d{l — 9). 

It remains to prove Proposition 4.4. 

Proof of Proposition 4-4- Using the notation of Definitions 4.2 and 4.3 we can 
find a finite number of points xi^n, ■ ■ ■, XL„,n G such that 

- xy„| > e^” for all 1 ^ f j < L„, (4.23) 

where in fact := L„(0) := |n„(0)|. In particular, (4.18) ensures that 

a^ndii-e) ^ for all n. (4.24) 

It might help to recall that a G (0,1) is independent of the parameters {i ,j , n) 
of this discussion. 

Even though F might not support a particularly-natural measure, the fol¬ 
lowing defines a very natural locally finite measure p, on E: 

OO Ln 

(4.25) 

n—M j=l 

for all F C R'^. 

Because ^(5„) = F„, the total ^-mass of the nth shell satisfies 

aQr^dii-e) ^ ^ n^M. (4.26) 
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Since /i is a measure on E, we plan to use the measure /x in Lemma 2.5 in order 
to find a lower bound for Dimjj(i?). With this aim in mind, we next establish 
an upper bound for ijl{Q) for every uprightbox Q whose sidelength is at least 
one, with the sole provision that Q C Sn for some n ^ M that is sufficiently 
large. 

Let us consider an arbitrary upright box Q{z ,r) of sidelength r ^ 1 such that 
Q{z,r) lies entirely in Sn for some large enough integer n ^ M. Elementary 
properties of the Euclidean space imply that there exists a positive integer 
k ^ 2‘^ together with k points zi,... ,Zk from the collection ..., XLn,n} 

such that 

k 

Q{z ,r) C [j Q{zj ,r). (4.27) 

i=i 


Therefore, 


k 

fi{Q{z , r)) ^ ^ KQ{zj ,r)) = k ^ 2"*. 
i=i 


This shows that for all p > 0, 


(4.28) 


Kp := sup 


^^{Q) 


Q G B, Q C Snj side(Q) G [1 ,r], r ^ 1 ^ ^ 2 


(4.29) 


On the other hand, if r ^ e®", then an upright box Q{z ,r) can contain at most 
(1 + re“®")‘^-many points in n„(0). Now, 

(1 + re"®”)'^ < 2‘^r^e-‘^^”, (4.30) 

because r > e^", and side(Q( 2 ,r)) < e” because Q{z ,r) C 5„. Therefore, it 
follows from the construction of the measure /i that 

I r ■■ Q & >3, Q C Sn, side(Q) > e^”| 

I [side(g)]^ J 

^ 2d^-d0n ^d-p ^ 


as long as 0 < p < d. Now we compare (4.29) and (4.31), and set p := d{l — 9) in 
order to see that Kdd-g) < where Kp was defined in (4.29); see also (2.14). 
This bound and (4.26) together yield the following: For all n ^ M sufficiently 
large, 

^dii-e) (^) ^ 2-''e-"''(i-^)p(5„) > a2-\ (4.32) 

thanks to Lemma 2.5. It follows that limsup„_,.oo > a2~‘^ > 0, and 

hence 

Dimjj(E) >Dim^E^d(l-6i), (4.33) 

where Dinig denotes the lower Hausdorff dimension of Barlow and Taylor [3,4], 
as was recalled in (2.16). □ 
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Let us introduce a last piece of notation before we state and prove the general 
lower bound of this section [Theorem 4.7 below]. That lower bound will be a 
counterpart to Theorem 4.1. 


Definition 4.6. Let X denote the collection of all independent finite sequences 
of independent random variables. 

Then we have the following general lower bound statement. 

Theorem 4.7. Suppose there exists b G (0, oo) such that C{h) < oo. Suppose 
in addition that there exist 6 € ( 0 , 1 ) and an increasing nonrandom measurable 
funetion S' : IR —>■ R sueh that 


n ^ max max inf logP{|S(Xt ) — S(y,)| > 1} —>■ —oo, (4.34) 
as n ^ CO. Then, 

X / d \ 

a.s. Moreover, if j £ (0,(i) then 

Dim„ |t G T : |t| > exp(e), ^ (^) } ^ ^ 

Remark 4.8. Condition (4.34) is a coupling assumption, and states that if 
ti, ■ ■ ■ ,tm & T have large norms [say, are in Sn for a large n] and sufficiently far 
apart [say at least exp((5n) apart], then ... ,Xt^ are close—say within one 
unit—to an independent sequence with very high probability. At first glance 
this might seem to be a technical and complicated condition. We will see in the 
next few sections that (4.34) is in fact frequently easy to use, particularly in 
the context of stochastic PDEs. Condition (4.34) can be recast in terms of the 
“correlation length” of the process A; see Conus et al [13] for details. 


If Dimjj G > 0 then in particular G is unbounded. In this way we see that 
(4.36) implies that 

a.s. for all 7 G (0, d). Let 7 f d to deduce (4.35) from (4.36). Thus, we need to 
derive only (4.36). 


Proof of Theorem 4- Since At ^ a if and only if S'(At) ^ S(a), we can replace 
the random field {At}tgT by the random field {S(At)}tgT throughout the entire 
statement of the theorem in order to see that for the remainder of the proof we 
can—and will—assume without incurring any loss in generality that 


S{x) := X for all a; G R. 


(4.38) 


In other words, the function S merely plays the role of a change of “scale.” 
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We plan to prove that the random set —dehned earlier in (4.14)—a.s. is 
6 *-thick for every 6 G (y/c?,!) and 7 G (0,d). Owing to Proposition 4.4, this 
proves that 

Dimjj A£ ^ (i(l — 0) a.s. for all 0 G 7 G (0 > (4.39) 

and (4.36) follows. In light of the paragraph that followed the statement of 
Theorem 4.7, this endeavor completes the proof of Theorem 4.7. Henceforth, 
we choose and fix two arbitrary nnmbers j G (0,d) and 0 G (j/d, 1). We also 
hold fixed an arbitrary [small] 

O<d<0-^. (4.40) 

Now we carry out a multi-scale argument. Recall the definition of n„(0). 
Because of that definition, for all sufficiently large integers n ^ 1, we can find 
points xi^n, ■ ■ ■ ,XL„,n m such that 

Qixi,n , e^”) n Q(xy„ , e^”) = 0 when 1 < i j ^ L„, (4-41) 

and 

aend(i-0) (4.42) 

where a G (0,1) depends neither on n nor on the pair {i ,j). For future purposes, 
we would like to emphasize that as part of the construction of these points we 
have also the following: 

lim — — = (i(l — 0). (4.43) 

n—¥oc> Ti 

For all 1 ^ i ^ we can find points ■ ■ ■ 5 'Z£„(i),n (0 Q{xi,n ,e^") 

such that 

\zk,n{'i) - Zl,n{i)\ > (4-44) 

whenever \ ^ k ^ I ^ tn(i), and 

be^dio-s) ^ ^ ^-i^ud{0-s)^ (4 45 ) 

where b G (0,1) depends neither on n nor on the triple {i ,k ,1). In fact, a and 
b depend only on 0 and S, which are held fixed throughout this discussion. 

According to (4.34) for all AT > d and 1 < i we can find an independent 
sequence Yi,, Ym of random variables such that 

max ~^j\ > 1} ^ ATe”^" for all n ^ K. (4.46) 


The particular construction of Yi,..., Ym might—or might not—depend on K 
and i; it does not matter. The upshot is the following: Since 


P 


/ { 7 

\teQ(lZ,e^-) (logf)'/^ " \Cib) 


!/&' 


^ P < max X,. (i\ < 


( (n + l) 7 y/^| 

I C{h) ) /’ 


(4.47) 
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two back-to-back applications of (4.46) show us that 

r w. / xi/o 


p. 




[teQC 


( 1 

(™neO") (logt)!/'' ^ \C{b) 




1/b^ 


^ Ke' 


i=i 
->(i) 




Ke 


-Kn 


(4.48) 

(4.49) 


i— 1 \ \ y / 

uniformly for all 0 ^ i ^ L„ and n ^ K. Given an arbitrary e G (0,1), we can 
find Kq > K such that 


< 1 _ p-(7-l-£)" 


1 /b^ 


^ 1 - e-(^+^)”, (4.50) 

uniformly for all 0 ^ ^ L„, 0 ^ j ^ ^n(*), and n ^ Kq. This follows because 

C{b) < oo. 

Because of (4.40), we can—and will—in fact choose e small enough so that 

0 < e < (4.51) 


Therefore, there exists Ki > Kq such that 

J 


I I 

(teQG^n.e®'*) (logt)i/^ ^ VC'(&)4 


< ife-^” -f (l - e-('^+®)” -f Ke 


-Kn 


(4.52) 

,G(i) 


< Ke -I- exp 


/ 

{-^„(,)e-(7+2^)"}. 


uniformly for all 0 ^ ^ L„ and n ^ Ki. [The preceding hinges on the fact 

that K > d > 7 -I- £.] We may deduce from (4.45) that 

'.i/f-'l 

(4.53) 


P 


( ^ ^ { 1 
(logt)i/^ - {C{b) 


^ Ke-^^ + exp i^-be‘^d 0 -dS-^- 2 e)n^ ^ 

Ln and n ^ Ki. Thanks to (4.51) and the fact that 
ee 14.43')—it follows that 


uniformly for all 0 ^ z ^ Lri and n ^ Ki. Thanks to (4.51) and tl 
K > d> n~^ logL„ = d(l — 9) + o(l)—see (4.43)—it follows that 

EEpI »p 7^4Tir«(^)‘'‘l<»- 

a—li—0 


[teQG^n.eO") (logt)i/^ ^ VC'(&)4 


(4.54) 
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Therefore, the Borel-Cantelli lemma ensures that the following holds for all w 
in the probability space that lie outside of a single set of P-measure zero: For 
all but a finite number of integers n ^ 1 , 


sup 


Xtjcj) 

) {\ogty/b 


> 



for all 0 ^ i ^ 


(4.55) 


Recall the random set that was defined earlier in (4.14). We can state the 
preceding display in another way; namely, that for all but a finite number of 
integers n ^ 1 , 


Ai n Q{xi^n, e^") 0 for all 0 ^ i ^ L„. (4.56) 

This proves that is 0-thick a.s. Proposition 4.4 then shows that 

Dimjj > c?(l — 9) a.s. (4-57) 

Since 0 G (y/d, 1) were arbitrary, we let 0 I y/d to complete the proof. □ 

Let us close this section by recalling a well-known general asymptotic evalua¬ 
tion of the tail of the distribution of the supremum of a stationary Gaussian pro¬ 
cess with a nice covariance function. The result is originally due to Pickands [36, 
Lemma 2.5], with extra conditions that were removed subsequently by Qualls 
and Watanabe [37, Theorem 2.1]. Albin and Choi [2] contain a novel elementary 
proof, together with an indepth discussion of the literature of the subject. 

Lemma 4.9 (Pickands). Let {? 7 (t)}tyo denote a continuous stationary Gaussian 
process with JHijit) = 0 and Yar r]{t) = 1 for all t ^ 0. Suppose that there exist 
constants k G (0, oo) and a € (0,2] such that 

Cov[? 7 (t), 7 y( 0 )] = 1 — -I-o(t“) ast—(4.58) 


Then, 

P-f sup ? 7 (s) > xl = /2 ^ ^ (4.59) 

I se[o.i] J (27r)V2 

where Ha G (0, oo) is a numerical constant that depends only on a. 

Remark 4.10. The cited literature also contains the assumption that there 
exists h G (0 , oo) such that inf^gjo^/ij Cov[r]{t), 77 ( 0 )] > 0. We have omitted it as 
it is subsumed by the assumed behavior of Cov[? 7 (t) , 77 ( 0 )] near t = 0. 

Remark 4.11. The Pickands constant Ha is itself rather interesting. We follow 
Pickands [36] and let $ := {4>(t)}t^o denote a centered Gaussian process with 

Cov[$(s), $(t)] = s“ -p - Jt - sl“. (4.60) 
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The process $ is a fractional Brownian motion with parameter a/2 when a G 
(0,2), and ^{t) = for a standard normal random variable ( when a = 2. 
Then, is the following nontrivial limit [which is known to exist as well]: 


Ha = lim — E 
T->oo T 


sup e 

tG[0,T] 




(4.61) 


It is known that Hi = 1 and H 2 = tt Other values of Ha are not known. 
See Harper [22] and the references therein for recent estimates. 


5 Peaks of the Ornstein-Uhlenbeck process 

For a first, and perhaps simplest, example of the general theory of §4 let us 
continue to write B for a standard Brownian motion on R, and define 

U{t) := e“*/2^(e*) for all t ^ 0. (5.1) 


Then U := {U{t)}t^o is an Ornstein-Uhlenbeck process; that is, f7 is a centered 
Gaussian diffusion with E[{7(t)f7(s)] = exp(—|t — s|/2) for all s,t ^ 0. Thanks 
to the law of the iterated logarithm for Brownian motion, 


lim sup 

t—¥oo 


Ujt) 

(21ogt)i/2 


= 1 


a.s. 


Here we consider the exceedance times of U, defined as follows: 

^ 7 } (7>0). 


(5.2) 


(5.3) 


This notation is consistent with the notation in (1.7) and (1.8) with g{x) := 
(2 log+ x)^/^. 

Because Vuil) = ^ogTsil )—where Vb{i) was dehned in (3.1)—and the 
natural logarithm is strictly monotone, we see that the random sets Vu (7) and 
Vb(i) are bounded and unbounded together. In particular. Proposition 3.4 
implies that Vuil) is unbounded [a.s.] if and only if 7 < 1; this fact follows also 
from the integral test of Motoo [29]. Therefore, the following theorem computes 
the macroscopic Hausdorff dimension of Vuil) in all cases of interest. 

Theorem 5.1. Part (1.20) of Theorem I .4 holds. Namely, 

Dimjj Vuil) = 1 — 7 ^ a.s. for all 76 (0,1]. (5-4) 

Remark 5.2. We can compare Theorems 3.2 and 5.1 in order to see that 
Dimjj logPB(l) = 0 a.s. whereas Dim^ Pb( 1) = 1 a.s. Equivalently, 

Dimjj exp(Rc/(l)) = 1 ^ 0 = Dimjj R[/(l) a.s. (5.5) 
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The elegant theory of Weber [40, Theorem 6.1] implies the following closely- 
related result: With probability one, 

lim - log |{0 < j < e" : Vu{l) n [j, j + l)^0}\ = l- 7 ^, (5.6) 

n—>-oo 77, 

where j • • • j denotes cardinality here. In other words, 

Dim„P; 7 ( 7 ) = 1 - 7 ^ a.s., (5.7) 

where Dim,^, denotes macroscopic Minkowski dimension. In the notation of Bar- 
low and Taylor [3,4], Dim,^^ E is the common value of dimyM E and dimLM E, 
when the two are equal. Since Dim^ E ^ Dim^, E for all E C [3, Lemma 
3.1(i)], Weber’s theorem implies half of Theorem 5.1; i.e., Dim^ Vuil) ^1 — 7 ^ 
a.s. As part of proving the converse inequality, we plan to use Theorems 4.1 
and 4.7 in order to obtain both inequalities at the same time. 

Proof of Theorem 5.1. We apply Theorems 4.1 and 4.7 with b = 2, Xt := U{f) 
for all < ^ 0, and S{x) := x for all x S R. Since {U{f)}t^Q is stationary, 
we can see from an elementary bound on the tails of the Gaussian law that 
c(2) = C(2) = 2. 

Because 

Cov([/ (t) ,U{0)) = e“*/^ as s —>■ t, (5.8) 

Pickands’s lemma [Lemma 4.9] implies the maximal inequality (4.5), and our 
Theorem 4.1 then shows that Dimjj Vuil) ^1 — 7 ^ a.s. for all 7 G (0,1). 

In order to prove the reverse inequality, let us note that if < • • • < tm are 
in Sn and ti+i — U ^ expjdn} for all I ^ t ^ m, then we may set 

Ti := [i?(e**) — B(e**“^)] (1 ^ z < m), (5.9) 

with to ■= It is easy to see that the l^’s are independent and 

max E(\U{U)-Yif) = ^expl-e*^”}. (5.10) 

Consequently, a standard bound on the tails of Gaussian laws implies that 

max P{|[/(ti) — Til > 1} < 2exp (—i exp {e*^"}) . (5-11) 

Of course, this proves that 

lim — max logP{|{7(ti) — lij > 1} = — 00 , (5.12) 

n^oo Tl 

with room to spare. Hence, Condition (4.34) is verified since the Yfs are in¬ 
dependent. We can apply Theorem 4.7—with 7 replaced by 7 ^—in order to 
deduce that ^1 — 7 ^ a.s. for all 7 G (0,1). This completes the 

proof. □ 
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6 Peaks of the linear heat equation 

Now we move on to examples that are perhaps more interesting. Consider the 
linear stochastic heat equation, 

Zt{x) = \Z^{x) + ^t{x) (x G R, t > 0), (6.1) 

subject to Z[j{x) := 0 [say], where ^ denotes space-time white noise. That, ^ is 
a totally-scattered centered Gaussian noise with 

Cov(^t(x) ,^s(y)) = (5o(s - t)(5o(a; - 2 /) for s,t ^ 0 and x,y G R. (6.2) 

It is well known—see Walsh [39, Chapter 3] for example—that there exists a 
unique integral solution to the stochastic PDE (6.1), and that solution has the 
following representation in terms of a Wiener integral: 

Zt{x) = f pt-s{y - x)^{dsdy) (t>0, ,xGR), (6.3) 

d(0,i)xR 

where the function {s,t,x,y) :—>■ pt-s{x — 2 /)ll(o,oo)(^ ~ s) denotes the funda¬ 
mental solution to the heat operator, 

o.(0,oo)xK. (6.4) 

That is, 

- 7 ^' 

for every t > 0 and x G R. 

It is also well known—see Walsh [ibid.] —that the random field Z has a 
modification that is continuous in (t, x); from now on we always use that version 
of the random field Z in order to avoid measurability issues. 

We are interested in the structure of the tall peaks of the random field 
X :—>■ Zt{x), where t > 0 is held fixed. With this aim in mind, consider the 
random set 

P2.(7) := |x > exp(e) : S 7 (;) |. (6.6) 

where t, 7 > 0 are fixed. 

Theorem 6.1. Every T’zti'l) is almost surely unbounded ifj ^ 1; else, (fy > 1 
then Vzt ( 7 ) is almost surely bounded. Furthermore, 

Dimjj Pze ( 7 ) = 1 - 7^ a-Sn (6.7) 

for all t > 0 and 7 G (0,1]. 
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A moment’s thought shows that Theorem 6.1 is an equivalent reformulation 
of Theorem 1.3 of Introduction. From a technical point of view, this particular 
formulation is more natural to state as well as prove. 

As we will see very soon, the fact that 7 = 1 is critical for the unboundedness 
of Vzt ( 7 ) is a fairly routine consequence of well-known theorems about the 
growth of a Gaussian process [37]. The main assertion of Theorem 6.1 is the one 
about the macroscopic Hausdorff dimension of Vztil)- Still, let us mention also 
the following immediate consequence of the first [more or less routine] portion 
of Theorem 6.1: 


lim sup 

X—>-oo 


Zt{x) 

(2 logxy/'^ 



a.s., 


( 6 . 8 ) 


for every nonrandom t > 0. A “steady state” version of this fact appears earlier, 
for example, in Collela and Lanford [10, Theorem 1.1(c)]. The following lemma 
puts things in the general framework of Gaussian analysis. 


Lemma 6.2. Fix some t > 0. Then, {Zt{x)'\x^Wi is a stationary Gaussian pro¬ 
cess with EZi(O) = 0, VarZt(O) = (I/tt)^/^, and Corr[Zt(x), Zt(0)] = 0(ja;j““) 
as jxj —>■ 00 for every a > 0. Finally, 


1 /7r\i/2 

Corr[Zt(a;), Zt(0)] = 1 - - Jx] -I- odxj) as Ja;]-)> 0. (6.9) 

Proof. Clearly, x 1 —>■ Zt{x) is a mean-zero Gaussian process with 


Cov [Zt{x), Zt{x')] = f p2s{x — x')ds for all X, G R. (6.10) 

Jo 

It follows from this formula that Zt{») is stationary as well, and has variance 
VarZt(0)= / p2s(0)ds= / (47rs)“^/^ ds = . (6.11) 


Jo 


Jo 


Furthermore, the preceding display shows also that 

Cov[Zt{x), Zt{0)] = f p 2 six)ds (6.12) 

Jo 

is bounded above by a Hnite constant C{t) times and hence goes to 

zero faster than any negative power of jxj, as jxj —>■ 00 . Finally, we note that if 
X > 0, then 


Var(Zt(0))-Cov[Zt(x),Zt(0)] 


b2s(0) -P2six)] ds 

fit/. " 


4v^ JO 

X 

4-\/^ Jo 


(6.13) 

j r- 1/2 (1 _ dr 

r-1/2 (l-e-^/”) dr-hO(x 2 ), 
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as X 4- 0. It follows readily from this and symmetry that 


Corr[Zi(x), Zt(0)] = 1 — c|x| + O(x^) as x —>■ 0, (6.14) 

with 

r-i/2 (^1 - dr. (6.15) 

A change of variables shows that 

poo poo 

/ r-i/2(l-e-i/’')d?-= / s-3/2(l-e-")ds. (6.16) 

Jo Jo 

Write 1 — e“® = e“^ dy and apply the Tonelli theorem in order to see that 

c = 1 / 7 ( 741 . □ 

We are ready to establish Theorem 6.1. 

Proof of Theorem 6.1. Throughout the proof, we hold fixed an arbitrary t > 0. 

Lemma 6.2 verifies all of the conditions of Theorem 1.1 of Qualls and Watan- 
abe [37], and hence it follows from that result that Vztil) is a-s. bounded if 
7 > 1 and a.s. unbounded if 7 ^ 1. In particular, we obtain ( 6 . 8 ) immediately. 
Furthermore, we can see—using the notation of §4—that 

6 = 2 and c( 6 ) = C(h) = • (6-17) 

Thanks to Lemma 6.2, the assumptions of Lemma 2.5 of Pickands [36] are 
met. Lemma 4.9 [Pickands’ theorem] implies the maximal inequality (4.5); 
therefore, we may apply Theorem 4.1—with 7 replaced by 7 ^—in order to con¬ 
clude that 

Dimjj VzAl) < 1 - 7 ^ (6-18) 

a.s. for all 76 (0,1]. We plan to prove a matching lower bound by appealing 
to Theorem 4.7 with S{x) := x for all x G R. Therefore, it remains to verify 
the coupling assumption (4.34), which we do next. 

For every i? > 0 we may define a space-time Gaussian random field as 
follows: For all x G R, 



■= f Pt-siy - x)^{dsdy). (6.19) 

Ji0,t)x[x-(Btp/2, a;+(Bt)i/2] 

It is intuitively clear that that Z « Z^^^ when i? ^ 1. Next we claim the 
following quantitative improvement of this remark: For all t, B, X > 0, 


sup P 





( 6 . 20 ) 
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Indeed, because Ps{z) < Ps(0) = (27rs) for all s > 0 and 2 G R, the Wiener 
isometry yields 


E 


Ztix) - Zj^\x) 


= [ ds [ dz [Psiz)]^ 

Jo i|z|>(Bt)i/2 


< 


ds 

\^2tts 


P 






( 6 . 21 ) 


where X has a standard normal distribution. If s G (0 , t), then we combine the 
elementary bound, 


P 




( 6 . 22 ) 


with a standard bound on the tails of the standard normal distribution in order 
to see 


Var (^Zt{x) - Z^^\x)'j < 



(6.23) 


The claim (6.20) follows readily from this and another appeal to the tails of the 
Gaussian laws. 

We use (6.20) in order to prove (4.34) using the following. 


Observation 1. If xi,X2, ■■■ ,Xm G IR satisfy \xi — Xj\ > when 1 < 

i ^ j ^ m, then the random variables zj:^\xi ),..., zj:^\xm,) are independent. 


Choose and fix some d G (0,1). If e” ^ xi < • • • < Xm < are m 
arbitrary points in such that Xi+i — Xi ^ exp{5n}, then we set Yj := z[^'^ (xj) 
for all I ^ j ^ m. Thanks to Observation 1, Yi,, Ym are independent random 
variables as long as n is large enough to ensure that 2\/nt < exp{(In}. And (6.20) 
ensures that 

max P{|Zt(xi) - y*| > 1} ^ 2exp (6.24) 

\ Z \ ot J 

In particular, 

lim — max logP{|Zt(a;i) — > 1} = —oo. (6.25) 

n^oo n l^i^m 


This implies (4.34) readily, and the lower bound that complements (6.18) follows 
from the conclusion of Theorem 4.7. □ 


7 Peaks of a nonlinear stochastic heat equation 

Let us now consider the following nonlinear stochastic partial differential equa¬ 
tion, 

ut{x) = \u'f{x) Y a{ut{x))ft{x) (x G R, t > 0), 

uq{x) = 1, 
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(7.1) 

(7.2) 



where ^ denotes space-time white noise, as before, and cr : ]R —>■ R is a Lipschitz 
continuous and non-random function with ct(0 ) = 0. 

It is well-known that the stochastic heat equation (7.1) has a unique solution; 
see Dalang [16], for instance. And that solution solves the following stochastic 
integral equation, interpretted in the sense of Walsh [39]: 


ut{x) = l+ pt-s{y - x)a{us{y))^{dsdy), (7.3) 

d(o,t)xR 

where pt{x) is the standard heat kernel on (0 ,oo) x R; see (6.5). 

It is known also that the solution to (7.1) is strictly positive for all t > 0; see 
Mueller [30] for a closely-related statement. The precise positivity assertion that 
is required here follows from the work of Mueller and Nualart [31]. Therefore, 
the tall peaks of a; i— Ut(x) and x i—>■ ht{x) match, where 

htix) ■■=logUtix). (7.4) 


The random field h is particularly well studied when a{z) = z for all z G R. 
In that case, a formal change of variables suggests that 

A*(x) = ih"(x) + i(h((x))2-6(a:), (7.5) 

subject to ho{x) = 0. This purely-formal “computation” is analogous to the 
classical Hopf-Cole solution to Burgers’ equation, and is due to Kardar et al [26]. 
The resulting ill-posed stochastic PDE (7.5) is the so-called “KPZ equation” of 
statistical mechanics. The recent solution theory of Hairer [21] gives a meaning 
to the analogous version of (7.5) where the x variable lives in [0,1] [together 
with suitable boundary conditions]. As far as we know, the original problem on 
R has not yet been given a rigorous meaning. 

In this section, we plan to study the set of points x > exp(e) at which the 
solution ut{x) exceeds certain high peaks. For the parabolic Anderson model— 
that is when a{z) = z for all z G R—Conus et al [12] have demonstrated that, 
for every t > 0, the tall peaks of a; >->■ ht{x) are of rough height (log|a;|)^/^ as 
I a: I —>■ oo. Specifically, they have proved that 


0 < limsup 

X—¥00 


htjx) 

(loga:)^/^ 


< OO 


a.s., 


(7.6) 


for all t > 0. Conus et al [ibid.] have also proved that the function (loga;)^/^ fails 
to correctly gauge the height of the tall peaks of ht {x) for general nonlinearities 
a. 

We will prove among other things that (7.6) holds when |(t(z)/z| is bounded 
uniformly from below by a positive constant. The mentioned boundedness con¬ 
dition is known to be an intermittency condition for the system (7.1) [18]. 

In order to describe our results in greater details let us define 


C := inf |ct(z)/z|, To- := sup |cr(z)/z|. 
zeR\{o} zeR\{o} 


(7.7) 
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Because a is Lipschitz continuous we always have 0 ^ £a- ^ La- ^ oo. We will 
be assuming that 



0 < .^cr ^ Lcr < 00 . 

(7.8) 

We call (7.8) an 

“intermittency condition” because it is the only known condition 

under which the solution to (7.1) is known to be intermittent 

in the sense that 


, Hk) . . , . . r,. 3 

k !->■ —j— IS strictly increasing on [ 2 , 00 ), 

fC 

(7.9) 

where either 

A(fc) := limsupt“^ supE ( Mt(x)/') , 

(7.10) 


t^oo a:eR 


or 

\{k) := liminf t~^ inf E ( U((a;)/') , 

(7.11) 


describe respectively the top and bottom fcth moment Lyapunov exponents of 
the solution; see Foondun and Khoshnevisan [18, Theorem 2.7].® 

Now define for all f, 7 > 0, 


VhAl) ■■= ja; ^ exp(e) : • (7.12) 

We will use the general theory of §4 in order to prove the following, which is 
the main result of this section. It might help to recall yet again our earlier 
convention that when we state that Dim^ E < 0 we mean that E is bounded. 

Theorem 7.1. Under (7.8), the following holds with prohahility one: 


(I) ' ^ 'TJS' « (I) ' (7.13) 


for every t ,7 > 0. Moreover, 

1 — ^ Dimjj Phtil) ^ 1 — (97^^^ 

where Phtil) "was defined in (7.12), 

4^2 . ^ 4V2 

“:=W' 


(7.14) 


(7.15) 


When a{z) = z for all z £ H, the stochastic PDF (7.1) simplifies to the 
following, which is known as a parabolic Anderson model and/or diffusion in 
random white-noise potential: 


iitix) = ^u'fix) -U Utix)ftix) (x G R, t > 0); 
uoix) = 1 . 


(7.16) 


In this case. Theorem 7.1 yields the following exact formula, which is an equiv¬ 
alent but perhaps more explicit way to state Theorem 1.2. 

®In general, the Lyapunov exponents, as were describe in the introduction, do not exists. 


33 



Corollary 7.2. The solution u to (7.16) satisfies the following: For all 7 , 1 > 0, 


Dinijj ^ exp(e) : Ut{x) > ^ (log®) ^ | = 1 _ 


4 ^ 3/2 

3 ' 


a.s. (7.17) 


Let us mention a rather general corollary of Theorem 7.1 as well. 

Corollary 7.3. Under (7.8), the tall peaks of the solution to the SPDE (7.1) 
are almost surely multifractal. 

Proof. Recall (7.15), let 71 := (2q:)“^/^, and then define 


7*+i := 



(7.18) 


iteratively for alH ^ 1. Clearly, 

1 - < 1 - = 1 - for all i > 1. (7.19) 

In addition, 0 < 1 — < 1 — < 1 for alH > 1 . 

Consider the collection of tall peaks Thtil) of order 7 G (0,1) that was 
defined in (7.12). Theorem 7.1 implies that 


Dimjj VhAli) < DiniH PhAli+i) for aU i > 1, a.s. (7.20) 

Definition 1.1 then shows that the tall peaks of ht —whence also Ut —are a.s. 
multifractal. □ 


We begin the proof of Theorem 7.1 with a basic tail probability estimate. 
Proposition 7.4. For any t > 0, we have 


a 


liminfz inf logP{ht(a:) > 2 } > —- 

z—^oo tcGR y/t 

limsup sup logP {ht(a;) ^ z} ^ 

2—>-oo tcGR 


(7.21) 


where the constants ct and /3 were defined in (7.15). 

Proof. Let and respectively denote the solutions to (7.1) with a{z) := 
iaZ and a(z) := L„z. The moment comparison principle of Joseph et al [24, 
Theorem 2.6] tells us that because of the condition (7.8), 


E 



^ E ([Mt(x)]'=) <E 



(7.22) 


for all real numbers t > 0, a; G R, and fc ^ 2. We now use the Hrst part of 
Theorem 5.5 of Chen [ 8 ] and the Gartner-Ellis theorem, for example in the 
form of Chapter 1 of the recent book by Chen [9], in order to complete the 
proof. □ 
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Armed with Proposition 7.4 we can prove half of Theorem 7.1 quickly. The 
second, harder, half will require work that will be developed afterward. 

Proof of Theorem 7.1: Dimension upper bound. Our immediate goal is to es¬ 
tablish the dimension upper bound; that is, we wish to demonstrate the follow¬ 
ing: 

Dimjj Vht ^1 — 7 a.s. for all 7 > 0. (7.23) 

We claim that for every 7 G (0,1), 

supp/ sup ht{y) > f^logs^ 1 as s —>■ 00 . (7.24) 

zeR Iy6[z,2+1] VP / I 


If this were so, then it would show that Condition (4.5) holds with 5 = 3/2 
and Xx := ht{x), and Theorem 4.1 then implies (7.23). Thus, the dimension 
upper bound of Theorem 7.1 follows once we prove (7.24). According to (7.4), 
it remains to prove that for every 7 G ( 0 , 1 ), 


sup P < sup ut (y) > exp 
I yG[2,2+l] 



S-7+o(l)^ 


(7.25) 


as s —>■ 00 . 

Recall that, as a corollary to Proposition 7.4, we have the following slightly 
weaker variation on the desired estimate (7.25): If g : —>■ P+ is a nonrandom 

function that satisfies lims_>oo g{s) = 0 , then 


supP < ut{y) > exp 
yeR I 


f i-g{s) 

\ /3 



2/31 


< S-7+o(l)^ 


(7.26) 


as s —>■ 00 . 

In order to derive (7.25) from (7.26) we apply a chaining argument. With 
this in mind, let us first observe the following, which is a quantitative form of 
the Kolmogorov continuity theorem [28, Theorem C. 6 , p. 107]: There exists a 
Hnite constant r = r(t) > I such that for all real numbers fc ^ 2 


sup E 


sup 

x,x' ^[w 
x^x' 


\ut{x) 


Ut{x')\^'^ 

x'\k/2 


< re 


rk^t 


(7.27) 


This and the Chebyshev inequality together imply that uniformly for all real 
numbers r],s € ( 0 , 1 ), s ^ e, and k 2, 


sup P 


sup \ut{w) - Ut{y)\ > exp (vylogs)^/^ 

y^[w,w-\-£] 


< exp (^k^t — 2k{r] log s 


(7.28) 
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We apply the preceding bound with the following choices of e and k: 


/ 2 \ 1/2 

fc = A:(s):=f—j (r/logs)!/^. 

Let s*(? 7 ,t) := exp (?7“i[2Tt]^/2) to see that fc > 2 if and only if s ^ St{r],t). 
We apply (7.28) with these choices of £ and k in order to see that 


(7.29) 


supP< sup |■ut(^(;) - Mt(y)| > exp (? 7 logs)^/^ 

wGM. I y^[w,w-\-£] 


€ TS 


(7.30) 


uniformly for all rj G (0,1) and s ^ max{e, s*(? 7 , t)}. In particular, we might 
note that 


sup P < sup ut (y) > exp 


^ TS + sup P < Ut{w) > exp 
< s“')'+°(i) as s —>■ oo. 


logs 


/3 


2/3' 


logs 


2/3' 


— exp 


(?7 logs)^/^ 


(7.31) 


owing to (7.26). Every interval [z, z + 1] can be covered by at most e ^ + 1 
intervals of the form [w , m + e]. Therefore, the preceding implies that 


supP < sup ut{y) > exp 

I yG[ 2 , 2 +l] 


2 logs 


2/3' 


[£"1 + 1] s"^+°(l) = S"^+°(1^ 


as s —>■ oo. This proves (7.25), and hence the upper bound on the macroscopic 
Hausdorff dimension in Theorem 7.1. □ 


Now we begin to work toward establishing the lower bound in Theorem 7.1. 
In order to do that we will attempt to verify the coupling condition (4.34). A 
first attempt might be to follow the case of linear SPDEs/Gaussian processes. 
More concretely, we may try to follow the proof of the lower bound of Theorem 
6.1 and consider, for every B > 0, a space-time random field as follows: 
For alH > 0 and a; G R, 


.(B) 


(x) 


:= 1 + / pt-siy - x)a (u'f^y)) ^{dsdy). 


(7.32) 


It is not hard to apply a fixed-point agrument in order to prove that the random 
integral equation (7.32) has a unique solution u(^\ Moreover, that {t,x) !->■ 
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i 

u\ ' (x) has a continuous modification. [We will not prove any of this here since 
we will not need to.] 

The random field is the analogue of the random field that was 

defined earlier in (6.19), but we now interpret the stochastic integral in (7.32) 
in the sense of Walsh, whereas the one for can be understood in the sense 
of Wiener. 

The random field Z^^^ was introduced in the proof of Theorem 6.1 because 
Z^-^^ has the following two desireable properties: 

(i) Z''^^ Ri Z if B is large [see (6.20)]; and 

(ii) zj:^\xi),..., zj:^\xm) are independent if the Xj’s are sufficiently far 

apart from one another; for example, if \xi — Xi+i \ > [see Obser¬ 

vation 1.] 

By analogy, we might hope that: 

(iii) « M if i? is large; and 

(iv) u[^\xi ),..., u[^\xm) are independent if the Xi’s are sufficiently far apart 
from one another. 

If so, then we could use —in a similar way as we used Z^^'> —in order 
to verify Condition (4.34), thereby obtain a lower bound on the macroscopic 
dimension of the set of high peaks of u. 

As it turns out, (iii) continues to hold. However, (iv) is manifestly false; it 
is possible for example to show that the covariance of u). (xi) and u). {X 2 ) is 

strictly positive, for all xi,X 2 G H, no matter how far apart xi and X 2 are from 
one another. 

We remedy the situation by defining the following random fields instead: 
Choose and fix an integer H ^ 1, as before, and define 

u[^''^\x) := 1 for alH ^ 0 and a; G R. (7.33) 

Then we define random fields for every j ^ 1, iteratively, as follows: 


(R,m) 


(x) :=!+/ 

Ja 


(0.t)x[£c-(Bi)i/2, a:-|-(Bi)i/ 2 ] 


Pt-siy-x)a ^{dsdy), 


for every m ^ 1. The object of interest to us is the random field The 

following estimate shows that u « when B 1. 

Lemma 7.5 (A coupling lemma). There exists a finite constant K > 1 such 
that for all real numbers t > 0 and A > 1, and all integers B ^ 1, 


sup P 



u\^'^\x) >A|<Arexp 


(R + logA)3/" \ 

KVt J ■ 


(7.34) 
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Proof. According to Lemma 4.3 of Conus et al [12]—see also Lemma 10.10 of 
Khoshnevisan [28]—there exists a finite constant c such that 

supE( ut{x) — ^ cexp — il/c), (7.35) 

a:eR \ / 

uniformly for all real numbers fc ^ 2, A > 1, and t > 0. Therefore, Chebyshev’s 
inequality shows that 

P ^ Ut{x) — u[^'^\x) > a| < c inf exp — (log A + i3)A:) , (7.36) 

uniformly for all real numbers t > 0 and a; G R, and for all integers B ^ 1. This 
readily implies the lemma. □ 

Our next lemma is essentially due to Conus et al [12], and shows that the 
random process x i— u[^’^\x) decouples fairly rapidly. 

Lemma 7.6. Suppose t > 0 is a real number, B ^ 1 is an integer, and 
xi,...,Xm are points in IR such that \xi — Xj\ > 2B^^‘^y/t whenever i ^ j. 
Then, u[^'^\xi),... ,u[^'^\xm) are independent random variables. 

Proof. We include a proof for the sake of completeness. 

First, let us observe that if \xi — Xj\ > 2(51)^/^, whenever i ^ j, then 
t 6 (^’^^(a;i),... ,u\^’^\xm) are independent. This is because: (i) 

= 1+ [ pt-s{y - x)a{l)f{dsdy)-, (7.37) 

and (ii) If ipi,... jipm G i^(R+ x R) are nonrandom with disjoint support 
then the Wiener integrals f ip j df (1 ^ j ^ m) are independent [compute co- 
variances]. Next we apply induction, using the following induction hypothe¬ 
sis: Suppose that whenever \xi — Xj\ > 2f(i?t)^/^ for i ^ j, the random vari¬ 
ables u[^’^\xi),... ,u[^'^\xm) are independent. Then we wish to prove that 
if \xi — Xj\ > 2{£ + l)(Rt)^/^ [i ^ j], then ... ,u[^’^'^^\xm) are 

independent. This property follows readily from the properties of the Walsh 
stochastic integral; namely, that if < 1 >^,..., are independent predictable ran¬ 
dom fields and ip^,..., ip"^ G L^(R+ x R) are nonrandom with disjoint support, 
then the Walsh integrals f ipj<l>^ d^ are independent [1 ^ j ^ mj. This completes 
our induction argument, and prove the lemma. □ 

We can now verify the dimension lower bound of Theorem 7.1. 

Proof of Theorem 7.1: Dimension lower bound. Our proof of Theorem 7.1 will 
be complete once we demonstrate that 

Dimjj ^[ 7 / 0 :]^/^^ ^1 — 7 a.s. for all 7 > 0. (7.38) 
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In order to establish this fact we will appeal to Theorem 4.7; therefore, it remains 
to verify Condition (4.34). We will appeal, as we did in the proof of the upper 
bound, to the general theory of Section 4, using the identifications := ht{x), 
b := 3/2, and S{x) := exp(x). 

Let us choose and fix an integer n ^ 1 and a real number 6 G (0,1), and 
consider an arbitrary collection {xi}'^^ of points such that: (a) e” ^ xi < • • • < 
Xm < e"+^; and (b) Xj+i — Xi ^ exp((In). From now on, we set B := and 

y, := \ogu[^'^\xj) (1 ^ j < m). (7.39) 

According to Lemma 7.6, Li, • • • , Fm are independent as long as 2n^t < exp((5n); 
and Lemma 7.5 ensures that 

imaj^P {|5'(X^J - 'S'(y,)| > l} < AT exp ■ (7-40) 

[Recall that S{x) := x here.] Since the constant K does not depend on the choice 
of xi,... ,Xm G Sn, we have shown that Condition (4.34) holds. Theorem 4.7 
implies (7.38), and completes our proof of Theorem 7.1. □ 


8 d-dimensional diffusion in random potential 

We will conclude this paper by presenting examples of stochastic PDFs over 
R+ X IR,'^ where d ^ 1 is an arbitrary positive integer. In order to keep the 
ensuing theory at a reasonably-modest technical level, we focus only on linear 
stochastic partial differential equations of the following type: 

Ut{x) = ^{Aut){x) + Ut{x)T]t{x) for t > 0 ,a; G 
uo[x) = 1; 


where the Laplace operator A acts on the variable x G R'^, and rj := {??i(a:)}i^o,a:eR<i 
is a centered generalized Gaussian random field with covariance measure 

Cov[r]t{x),T]s{y)]=So{t-s)f{x-y) [s,t ^ 0, x,y G R'^j, (8.2) 

for a positive-definite bounded and continuous function / : R'^ —>■ R+. 

By comparison to the classical situation of Brownian heat baths, we can 
think of the solution to (8.1) as describing Brownian motion in a random envi¬ 
ronment. 

When d = 3, a variation of equations of the form (8.1) was introduced in 
cosmology in order to describe the large-scale structure of the universe. Here 
are some more details: It is believed that,® after we make a standard change of 

®See Albeverio et al [1] and its detailed references to the physics literature, in particular, 
to the pioneering work of Zel’dovich and his collaborators. 
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variables to remove mathematically-uninteresting physical constants, the veloc¬ 
ity held V of galaxy masses approximately solves the 3-dimensional stochastic 
PDE,'^ 

^vt{x) = ^{Avt){x) - (vtix))^ Vvt{x) + (8.3) 

for {t ,x) £ {0 , oo) X subject to the following: 

curl(dt(x)) = 0, vo{x) = -V-ipix); (8.4) 

for all t > 0 and a; G R^. The external held $ is a scalar held and believed to 
be random, and the initial held ^|) may or may not be random. Now we apply a 
formal Hopf-Cole transform and posit that 

vt{x) = -Vlog(/)t(x), (8.5) 

for a scalar held 4>. It is then easy to see that, if ip and $ were smooth, then p 
would solve 


(pt{x) = ^{A(pt){x) +(pt{x)^t{x) for t > 0,x G R^; 
cPo{x)=e^^^\ 


( 8 . 6 ) 


As far as we know, there is no general agreement on what the external held 
$ should be, though the simplest form of the “big-bang theory” might suggest 
that Ip = 6o, after a suitable relabeling of R^. Our stochastic PDE (8.1) is this 
equation in the particular case that $ = 77 and ip = 0. 


8.1 The main result 

It is a classical fact that our correlation function / is uniformly continuous and 
maximized at the origin; i.e., 

f{z) ^ /(O) for all 0 G R'^. (8.7) 

Therefore, in order to avoid trivialities we will always assume that 

/(O) > 0. (8.8) 

Indeed, if /(O) were zero, then / = 0 and hence 77 = 0. In that case, the solution 
to (8.1) is Ut{x) = 1, trivially. 

Let / denote the distributional Eourier transform of /. Because / is assumed 
to be positive dehnite, / is a positive distribution. That is, / is a tempered Borel 
measure on R'^ thanks to the Riesz representation theorem. Furthermore, the 
Parseval identity shows that 

/ f{<ix) = {2Trf [ {(p*if){z)f{z)dz^{2Trff{0)\\(p\\li^^ay (8.9) 

Jr-i Jr'I 

^Here, Vt; denotes the matrix of all first derivatives of the coordinate functions of v. 
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for all rapidly-decreasing test functions : IR,'^ —>■ R, where (p denotes the 
reflection of (p; that is, 

<p{x) := (p(—x) for all x G (8.10) 


We replace ip by pe —where {pe}e>o is an approximate identity built from func¬ 
tions in 5(R'^)—and let e 0 in order to conclude that / is in fact a finite Borel 
measure on R"^. In particular, this shows that 


[ fjdz) 

V 1 + lklP 


( 8 . 11 ) 


Thanks to the theory of Dalang [16], condition (8.11) implies the existence of a 
predictable mild solution to (8.1). Moreover, that solution is unique among all 
predictable solutions that satisfy 

sup sup E (|ut(x)|*) < oo for all T > 0 and k ^ 2. (8.12) 

ie[o.T] 


This solution to (8.1) can also be written in mild form as the a.s.-unique 
solution to the random integral equation. 


Ut{x) = {pt*uo){x)+ pt-s{x - y)us{y)r]{dsdy), (8.13) 

J(o,t)xn 

where the stochastic integral is understood in the sense of Walsh [39], and now 
p denotes the natural d-dimensional generalization to (6.5): 

e-ll^llV2i 

■= [i>0,xGR]. (8.14) 

Formally speaking, we can let d := 1, a{z) := z, and /(z) := So(z) to see 
that (8.1) is [in this case] one possible extension of (7.1) to higher dimensions. 
Although do is not a continuous and bounded function, it is an appropriate limit 
of such functions. As such, one can derive (7.1)—in this case—using a limiting 
procedure from the solution to (8.1) with a suitable approximate identity {/e}e>o 
in place of /. See Bertini and Cancrini [5] for the details. 

There is a good way to construct examples of positive-definite continuous 
and bounded functions / : R'^ —?> R+ as follows: 

f ■= h*h, (8.15) 

where h G L^(R‘^) is a nonnegative fixed function, and h denotes the reflection 
of h [see (8.10)]. Thanks to (8.7), we can see that 

/(O) = sup f{z) = \\h\\l 2 (^^^y (8.16) 

We will restrict attention to such correlation functions / only. In fact, we 
concentrate on a slightly-smaller class of correlation functions still. Namely, we 
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will consider only correlation functions / that satisfy (8.15) for a nonnegative 
h G L^( 1 R‘^) that satisfies the following: 

limsup^-^log ( [ [/i(2)]2dz 1 < 0. (8.17) 

n^oo logn \J\\^\\>n J 

It is possible to write a Feynman-Kac type representation of the solution to 
(8.1). That representation implies readily that Ut{x) > 0 a.s. for alH > 0 and 
X G With this remark in place, we have the following, which is the main 
result of this section. 

Theorem 8.1. Consider the SPDE (8.1) where the spatial correlation function 
f of the noise satisfies (8.8) and (8.15) for a nonnegative function h G L^(R‘^) 
that satisfies (8.17). Then, with probability one: 

limsup = ^/2tf{0)d, (8.18) 

IkiKoo (log||a:||)i /2 

and 

2 

Dim„ {x G R'" : ||x|| > e^ ut{x) ^ e^^i^s H-H | = d- (8.19) 

Recall that for all sets E, Dim^ E < 0 means that E is bounded. It follows 
easily from this convention that the lim sup law (8.18) is a consequence of (8.19). 
Conus et al [12] proved that the lim sup in (8.18) is strictly positive and finite 
a.s., and the evaluation of the lim sup is contained in the recent work of Chen 
[ 8 ]. We have included (8.18) merely to highlight the fact that the tall peaks of 
Ut{x) are of order exp{ 7 A/t log ||x|l} for a constant 7 > 0, and hence that (8.19) 
is indeed a multifractal description of the tall peaks of Ut{x). 

The remainder of this section contains the proof of Theorem 8.1. 

8.2 Proof of Theorem 8.1: Upper Bound 

The proof of Theorem 8.1 proceeds by verifying the conditions of the general 
theory of §4. Thus, the proof is divided naturally into two parts: Proof of 
the dimension upper bound; and a separate derivation of the dimension lower 
bound. 

The dimension upper bound will be obtained by verifying Condition (4.5). 
Our first lemma is essentially a specialization of the proof of Proposition 4.4 of 
Conus et al [14]. The only new observation is that since here / is continuous, 
the constants in Proposition 4.4 can be computed explicitly. 

Lemma 8.2. Eor all t > 0, 

uniformly for all x G R"^. 
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Proof. If TO ^ 1 is an integer and k G [to,to + 1), then Jensen’s inequality 
assures us that 


||ut(a;)||Lm(n) < ||■u^(a:)||ifc(^) ^ ||ui(a;)||L™+i(n)- (8.21) 

Therefore, it suffices to prove that the lemma holds where the limit is taken 
over all integers fc —>■ oo. 

Proposition 4.4 of Conus et al [14] includes the statement that 

E(|ut(a;)|'=) ^ (8.22) 

for all real t ^ 0, x G R'’*, and integers k ^ 2. 

We develop a corresponding lower bound by following the proof of Proposi¬ 
tion 4.4 of [14], but use the additional hypothesis that / is continuous. 

According to the Feynmen-Kac formula for the moments of the solution to 
(8.1)—see Conus [11] and Hu and Nualart [23]—the /cth moment of the solution 
to (8.1) has the following representation: 


E(|rq(x)|'=) =E 



(8.23) 


where are independent Brownian motions on R'^. 

Choose and fix some e > 0. There exists d(e) > 0 such that f{x) ^ /(O) — £ 
whenever ||x|| ^ d(£)/2. Now let us consider the event 


rip!— \ cj G ri I 


max sup |[A:j(s)|[(w) < S{e) 
is;js;fese[o.t] 


(8.24) 


The moment formula (8.23) and the continuity of / together imply that 

E (iMt(x)l'’’) ^ exp I (/(O) - £)t| P(He) (8.25) 

pi sup ||Xi(s)|| < 1 

[^G[0,t/|5(£)|q 

uniformly for alH ^ 0 and x G R"^, and all integers k ^ 2. This implies that 

liminf logE (|■ut(a;)|^) > (/(O) — £)t/2, (8.26) 

fc—>oo: keZ ' ' 

and readily yields the desired lower bound since £ > 0 were arbitrary. □ 

We can now invert moments, exactly as was done in the proof of Theorem 
5.2 of Chen [8], in order to deduce the following. We skip the proof since it 
really follows the proof of [8, Theorem 5.2] almost exactly. 


= exp I (f(0) - £)t| 
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Lemma 8.3. Fix t > 0 and a bounded domain D C Then uniformly for 
all X G R'^, 


lim z ^ logP {logwt(x) > z} 


lim z ^ log P 



sup ut{y) 

y<^Q{x,l) 


-{2fm-\ 




z 


(8.27) 


Lemma 8.3 verifies the conditions of Theorem 4.1, from which we can deduce 
half of Theorem 8.1 readily. 

Proof of Theorem 8.1: Dimension upper bound. Lemma 8.3 verihes Condition 
(4.5) of Theorem 4.1 with X{x) := \ogUt(x), z := ( 27 /( 0 )t logs)^/^ and b = 2. 
It follows from Theorem 4.1 that 

Dim, {x G R-' : ||x|| > e^ Ut{x) > eV‘iogil-llJ ^d- (8.28) 

almost surely for every t, 7 > 0. This completes the proof of the dimension 
upper bound in Theorem 8.1. □ 

Remark 8.4. Let us record also the fact that the preceding proof required only 
that /i > 0 is in L^(R‘^); the extra regularity condition (8.17) was not needed, 
and (8.28) holds without that extra condition. 


8.3 A coupling of the noise 

We will have need for a particular construction of 7 that can be found essentially 
in the paper by Conus et al [14]. Let f denote a space-time white noise on R+ x 
R*^; that is, ^ is a centered generalized Gaussian random field with covariance 
measure 

Cov[^t(x),^s(?/)] =(5o(t-s)-(5o(x-y) [s,t > 0, x,j/G R^']. (8.29) 

We can construct a cylindrical Brownian motion B, using as follows: For all 
(p G L^(R‘^) and t > 0, define 

Btip) ■= [ ip{y)f{dsdy) and Boip) := 0. (8.30) 

Then {Bt{(p)}t^o^^^L 2 (^d) is a centered Gaussian random field with 
Cov [Bt{(pi ), Bt{ip 2 )] = min(s , t) ■ {pi , 

for all s,t ^ 0 and (pi,p 2 G L^(R‘^). Thus, we see that {Bt}t^o is a cylin¬ 
drical Brownian motion on L^(R'^) [see Da Prato and Zabczyk [15, §4.3.1]. In 
particular, we can recover the space-time white noise f from B by noticing that 

nt/+l 
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where the derivative is understood in the sense of random linear functionals; see 
Chapter 2 of Gel’fand and Vilenkin [19] for this topic. 

We can subordinate a large family of Gaussian random fields to the cylindri¬ 
cal Brownian motion {Bt}t^o as follows: For all €. L^(]R‘^) we define a new 
centered Gaussian random field by setting 


(8.33) 

where we recall f/' denotes the reflection of f/' [see (8.10)]. 

The random mapping {(p , 1 ^) i— {(p) is and linear, and the covariance 
structure of each Gaussian process B^"^^ is dictated by 


Gov 


{p 2 ) =min(s,t)- , p 2 * tp * 


L2(]R<i) 

for every s,t^0 and pi,p 2 ,'P G L^(R'^). In particular, (8.15) yields 


(8.34) 


Gov 


{p2) =mm{s,t)-{pi,p2*f)L2(Rd), (8.35) 


for every s,t ^ 0 and pi,p 2 G This is another way to say that the weak 

derivative dtB[^\x) is a particular construction of the noise rjt{x) dt. Since we 
are interested only in the law of the solution u to (8.1), and that law is by 
construction a function of the law of ? 7 , we may—and will—change probability 
space if we have to in order to construct the noise rj on the new probability 
space as follows: 

rit{x) -.=—B[’^\x) (t ^ 0,a; G R'^). (8.36) 

Thus, we are justified in using the following notation to denote the Wiener 
integral dip. 


/r+ 


<^s{y)v{dsdy) := 


/R+xR^^ 


<S>siy)d,Bd^\y) dy, 


(8.37) 


for all nonrandom functions $ G L^(R+ x R'^). Thus, in particular, we may— 
and will—think of the solution to the stochastic PDF (8.1) as the unique solution 
to the following stochastic integral equation 


Utix) = l+[ pt-siy - x)usiy)dsB^’^\y)dy, (8.38) 

d(0.t)x]R‘^ 

for all t > 0 and x G R'^, where the stochastic integral is understood in the sense 
of Walsh [39] and written using the notation introduced earlier. 

For us, an advantage of this construction is that, in this way, not only have 
a construction of ijtix) for our fixed function h, but we have in fact produced 
a coupling of ip dtB[^\x) which is a linear map and agrees with r]t(x) when 
Ip = h. 

For purposes of comparison, let us mention that the stochastic differential 
(ds dy) of Gonus et al [14] is the same thing as our mixed random differential 
dsBi'^^y) dy. 
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8.4 Proof of Theorem 8.1: Lower Bound 

In order to prove the lower bound on the dimension in Theorem 8.1 we plan 
to verify Condition (4.34) using an approach that has the same flavor as the 
proof of the already-developed lower bound of Theorem 7.1. There are some 
nontrivial differences in the proofs, however. Most notably, since the spatial 
correlation of the noise of the present section is not in general 0 even when x 
and y are very far apart, we need to do something more. To combat this issue, 
we first approximate the function ft, by a compactly-supported function hp, and 
then follow the approach used in the proof of the lower bound of Theorem 7.1. 
Our proof will follow the ideas of Conus et al [14] loosely. 

Define, for all x G and /3 ^ 1, 


h^(x) := h(x)g^(x) where 0/3(x) := H 

j=i ^ ^ / 

Every function ft^ is in L^(R‘^) and has compact support. In addition, 
converges to ft pointwise as /3 —>■ oo. Therefore, we ought to be able to construct 
approximations of Ut(x) by first approximating the noise dtBl^\x) dx by the 
noise dtB[^^\x) dx of the previous subsection. 

In order to simplify notation, let us write for every x := (xi,..., Xd) G R'^, 
t > 0, and /3 ^ 0, 


Jt(x;/3) := (0 ,t)x xi - /SVt, xi + /SVt 


X • • • X 


Xd - pVt ,Xd + pVt 


We can follow the lead of Conus et al [14], and choose and fix /3 > 
consider the solution to the following stochastic integral equation: 


•'^^\x)-.= l+ [ pt-s{y - x)uf\y)dsB^s^^\y)dy. 
JXt{x-,l3) 


(8.40) 
0, and 


(8.41) 


The following was pointed out in [14, §5] without proof. 

Proposition 8.5. The stochastic integral equation (8.41) has a predictable so¬ 
lution that is unique among all such solutions that satisfy the following for all 
real numbers k 2: 


sup E 





for all P > 0 and t ^ 0. 


(8.42) 


We will need some of the ingredients of that proof, and the details are not 
included in [14]; therefore, let us hash out a few of the standard details. 

Keep P > 0 fixed, and define for alH ^ 0 and x G R"^, 

uP’°)(x) :=1. (8.43) 
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Then we define iteratively the random field as follows: For all integers 

m > 0, reals t > 0, and x € IR"^, 

:= 1 + / pt-s{y - x)uf''^\y)dsBf‘‘'\y)<iy. (8.44) 

JXt{x-,l3) 

The proof of Proposition 8.5 requires two a priori bounds. The first controls 
the moments of the Picard iterates. 

Lemma 8.6. Uniformly for all t ^ 0, x € k ^ 2, and integers m ^ 0, 

E (8.45) 

Proof. Throughout, we hold f3 > 0 fixed and define for all t > 0, m > 0, and 

X € R"^, 


Pt-s{y - x)uf''^\y)dsB^J^^\y)dy, (8.46) 

JXt{x-,0) 

if and when the stochastic integral is defined in the sense of Walsh [39]. As an 
a priori part of this proof we will derive moment bounds for Xf'^\x). 

Let 

fl^ := hg *hi3, (8.47) 

where hp was defined in (8.39). The starting point is that a suitable form of 

the Burkholder-Davis-Gundy inequality [28, Theorem B.l] implies that for all 
real numbers fc ^ 2, 

x\'^\x) ^ 4fc f T{s , y , z) dsdzdy 

./2 (8.48) 

sfdkfdsf dy [ dz ||T(s , j/, 2 ;)||fc/ 2 , 
do dR<i d]R<i 

where we have used Minkowski’s inequality in the last line, and 

T{s,y,z) := pt-siy - z)pt-siz - x) uf''^\y)uf''^\z) fyiy-z) 

(8.49) 

^ Pt-siy - z)pt-siz - x) uf^^\y)uf^^\z) fid), 

thanks to (8.7) and the elementary fact that hg ^ h, whence fp ^ /. In 
particular, by the Cauchy-Schwarz inequality, 

I|T(s, 2 /,z)|U /2 ^ 6^“® ATa.fc Pt-siy-z)pt-siz-x)fid), (8-50) 

for every a > 0, where for all space-time random fields dr := {'I's(2/)}s>o,y6]R‘^: 

A/'a.fc(4') := sup sup [e""'* ||^'s(2/)||fc] • (8.51) 
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We plug the preceding into (8.48) to see that for all a > 0, 


< 4fc " /(O) • e"“^ ds 

Ar„,fc (u(^'"*))]"/(0). 


2fc/(0)e 


(8.52) 


Multiply both sides by exp(—2at) and maximize over x G to see that 

Afc.,k ^ Afc.,k (8.53) 

Since di* h->- Afa,k{'^) is a bona fide norm on random fields, it follows from (8.44) 
and the triangle inequality that 

AT^.fc < 1 + (8.54) 

The preceding is valid for all a > 0. We now select a := 8fc/(0) in order to see 
that, for this special choice, 

■X'skf( 0 ),k ^ 1 + \X's,kf{ 0 ),k . (8.55) 

Since Na,kiu^^'^^) = 1 for all a > 0, induction implies that 

Ar8fc/(0),fc ^ 2. (8.56) 

This is another way to write the lemma. □ 

Next we state and prove the second a priori bound that is required for the 
proof of Proposition 8.5. 

Lemma 8.7. Uniformly for allt^0,xG R'^, fc > 2, and integers m ^ 0, 


E 


(/3,m+l) 


(x) - 


(0,m) 


(x) 


^ I ^ (8.57) 


Proof. The proof is similar to that of Lemma 8.6. Recall Xl^l’s from (8.46). 
By the Burkholder-Davis-Gundy inequality [28, Theorem B.l], 



2 

/ y{s ,y ,z)dsdzdy 


k 

JXt{x-,0) 


k/2 


if4k [ ds [ dy [ d^: ||3^(s , j/, 0 )||fc/ 2 , 
do dR<J d]R<i 


(8.58) 


where 

yis,y, z) := pt-s{y - z)pt-s{z - x)f{Q) 


(/3,m+l) 


{y) - 




(y) 


(P,m+1) 


iz)- 


,(d,rn) 


{z) 


(8.59) 
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Recall A/’a^fe '■= Afa,k from (8.51). Since 

Af8fe/(0).fe = Af8kfi0),k - x(-)) , (8.60) 

we can easily adapt the proof of Lemma 8.6 to see that 

J^8kfio),k < iAr8fc/(o),fc . (8.61) 

Because ul^’^\x) = 1, Lemma 8.6 implies that 

■^8kf{0),k < J^8kf(0),k +-^8fc/(0),fc ^ 3. (8.62) 

Therefore, Af8kf(o),k{u^^’™~^^'^ — < 3-2“’", which is another way to state 

the lemma. □ 

We are now ready to complete the proof of Proposition 8.5. 

Proof of Proposition 8.5. Recall Ma,k from (8.51). Lemmas 8.6 and 8.7 together 
guarantee the existence of a predictable random field such that 

■^ 8 kfi 0 ),k ^ 2 and Ji^ A/'8fc/(o),fc = 0. (8.63) 

The proof of Lemma 8.7 also implies that 

J5m^Afskm,k -X)=0, (8.64) 

where was defined in (8.46), and 

Xt{x) := [ pt-s{y - x)u[^\y)dsBf'^\y)Ay. (8.65) 

JXt{x-,P) 

These remarks together show that solves (8.41). Uniqueness is similar; in 
fact, the argument that led to (8.61) can be re-iterated in order to imply that if 
V were any other predictable solution to (8.41) that satisfies J^8kf(o),k(v) < oo, 
then 

X'8kf{0),k (v - < \J^8kf(0),k (r' - , (8.66) 

and hence M 8 kf( 0 ),k{v — =0. □ 

Now that we have justified the existence of a good solution to (8.41) we can 
establish that u « when j3 is large. 

Lemma 8.8 (A coupling lemma). There exists a finite eonstant K > 1 such 
that for all real numbers t > 0 and A > 1, and all integers /3 ^ 1, 

P {|ut(x) - > a} exp + ^ (8 67 ) 

uniformly over all x G 
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Proof. According to Lemma 5.3 of Conus et al [14], there exist finite constants 
c > 0 and b G (0,4) such that 

supE f ut{x) - ^ ^^ck^t-bkioep^ (g gg^ 

a:£]R \ / 

valid uniformly for all x € and all real numbers fc,/3 ^ 2 and t > 0. This 
bound and Chebyshev’s inequality together yield the following: Uniformly for 
all reals t > 0 and k,X,l3^ 2 , and x € 

P { Ut{x) - > a} < ^gCfc=^i-(logA+61og/3)fc^ (g gg) 

The preceding readily implies the lemma, after we optimize over k ^ 2. □ 

For every x, y G let us define 

H(a;,y) := miu jx; - 2 /il, (8.70) 

where we recall | • • • | denotes the £°°-norm on R'^. The following lemma is due 
to Conus, Joseph, Khoshnevisan, and Shiu. 

Lemma 8.9 (Conus et al [14, Lemma 5.4]). Suppose that t > 0 is a real number, 
13 1 is an integer, and x^^\ ..., are points in R*^ such that 

H > 2/3([log/3J + 1 ) ^1 +when 1 ^i ^ j ^ m. (8.71) 

Then, 

^(/3,Llog/3J + l) ^ ^ ^ ^(/3,Llog/3J+l) ^3^^2) 

are independent random variables. 

We can now verify the dimension lower bound of Theorem 8.1. 

Proof of Theorem 8.1: Dimension lower bound. Choose and fix a time variable 
t > 0. We will appeal to Theorem 4.7, specifically to the general theory of 
Section 4, using the identifications X(x) := logMt(x), b := 2, and ^(x) := 
exp(x). 

Let us fix a real number 6 € (0,1) and consider an arbitrary collection 
of points that satisfy the following: 

(a) x^*^ G C R*^ for all 1 < f ^ m'^-, and 

(b) Zl(x(*) ,x^^'>) ^ exp((5n) whenever 1 ^ f ^ ^ m'^. 

From now on, we set 

P := exp(n2/3) and Yj := y,-„ := log (^u^’LiogW+i) ^ (8 73 ) 
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for all 1 ^ i ^ m'^. We might observe that there exists N := Nt >0 such that 

2 ([log /3J + 1) /3 (l + Vi) < e-^” for all n^N. (8.74) 

Therefore, according to Lemma 8.9, Li, • • • , are independent whenever n ^ 
N ; and Lemma 8.8 ensures that 


max P 

l^j^m 


{|5(x(a:«)-5(r,) 



n4/3\ 

Kt J ■ 


(8.75) 


Since the constant K does not depend on the choice of ..., € Sn, 

we have shown that the coupling Condition (4.34) holds, with room to spare. 
Therefore, Theorem 4.7 implies that a.s.. 


Dimjj |a: G : ||x|| > e®, Ut{x) ^ | ^ d — 


2 /( 0 )’ 


(8.76) 


for all 7 > 0. In light of the already-proved upper bound this completes our 
proof of Theorem 8.1. □ 
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